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ABSTRACT 



A new adaptive control law for mechanical manipulators that maintains uniformly good 
performance over a wide range of motions and payloads is developed. This control strategy 
combines properties from both the Model Reference Adaptive Control and the Self Tuning 
Regulator Theory and serves to extend the Adaptive Model Following Control approach 
into using a nonlinear reference model. 

The design procedure is simple resulting in an overall system which is globally stable 
and offers itself to microcomputer implementation. The effectiveness of the approach is 
demonstrated on several computer simulations which compares its performances against 
some of the commonly known adaptive control techniques. 

Also presented is a comparison of the computation complexity of different methods used 
in deriving the dynamic equations of motion of a mechanical manipulator as well as a 
survey of various robot control methodologies available in the literature today. 
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I. ROBOTICS : AN OVERVIEW 



A. INTRODUCTION 

The 19S0's may easily be characterized as the robotics era. The last five or six years 
have experienced very strong social and economical demands for advanced automation in a 
fast expanding domain of applications ranging from the— well established car— making 
industry to unmanned underwater workstations [1]. Moreover, there is a widespread feeling 
that it is likely that robots, in the years ahead, will become crucial agents of industrial 
change, transforming production processes and affecting everyday lives [2], Confronted with 
these facts, we are led to wonder about the reasons behind experts attaching such heavy 
weight to robotics and about the characteristics that make the industrial robot such a 
powerful tool. 

An answer to these questions may be found by tracing the origins that tie robotics to 
automation. Roughly speaking, we can identify three types of automation [3]. 

1. Continuous process controls 

This type of automation is used in oil refineries. It employs mostly computers with 
no, or little, human intervention. This type is highly automated. 

2. Hard automation 

Hard automation uses mainly transfer conveyor methods to handle the high volume 
mass production. This type of automation is based on setting up specific assembly lines 
with special tools and gadgets. This implies that hard automation relies on hardware which 
cannot be easily changed, should a change in the design of a product be called for. 

3. Flexible automation 

Flexible automation handles low volume batch production. Because this type aims at 
overcoming the limitations of hard automation, it is also referred to as programmable 
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automation. This kind of "machine" is designed to be flexible and to be able to react to its 
environment in an adaptive fashion. This is clearly different from the conventional machine 
which can be used only for well-defined, specialized, and preappointed tasks. The 
mechanical manipulators used in industrial applications fall into this category. However, 
even though a considerable progress has been made in introducing robots into industrial 
situations, there is still more to learn both in overall concepts and practicalities before a 
robot having a performance comparable to humans can be built [4]. 

The next section will describe what is meant by a "mechanical manipulator", in 
general, and will outline some of its characteristics. 

B. A MECHANICAL MANIPULATOR: DEFINITION AND CHARACTERISTICS 
A robot is a computer-controlled mechanical device that can be programmed to 
automatically move objects through different configurations in space. Robots are normally 
constructed as series of coupled rigid links, which together constitute what is called a 
kinematic chain. There are two types of kinematic chains [5]: 

1. The linkage or the closed chain, where every link is connected to at least two other 
links in the chain. 

2. The manipulator or the open chain, where some of the links are connected to only one 
other link. The articulated portion of most industrial robots is an open kinematic chain 
with some form of end effector attached to the final link. 

A typical industrial robot is shown in Figure 1.1. The coupling of two adjacent links is 
referred to as a kinematic pair or joint. The most frequently encountered pairs in current 
industrial manipulators are the revolute or rotational joint and the prismatic or translation 
joint. These pairs are shown in Figure 1.2. The revolute and prismatic joints are single 
degree of freedom pairs. Any manipulator must have at least six degrees of freedom to 
enable it to achieve arbitrary real-world configurations. Thus, most industrial robots are 
constructed of exactly six links. The first three degrees of freedom are generally referred to 
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Figure 1 . 1 : A typical industrial mechanical manipulator 
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Revolute joint 



Figure 1 .2: A revolute and a prismatic pairs 
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as an arm subassembly. They are used primarily to position the wrist unit at the 
workpiece. The final three degrees of freedom are referred to as the wrist subassembly, 
subsequently employed to orient the tool according to the configuration of the object [6]. 
This research is mainly concerned with the arm subassembly. The orientation of the tool 
will not be considered. 

C. ROBOTICS APPLICATIONS 

Mechanical manipulators are widely used in manufacturing and assembly tasks such as 
material handling, spot and arc welding, parts assembly, paint spraying, loading and 
unloading numerically controlled machines, and in handling hazardous materials [7], 
Furthermore, it is now common belief that robot systems can be used in areas other than 
assembly tasks. Perhaps the most unusual application, to date, can be found in Australia, 
where a robot is used for sheep shearing [8]. Another application of robots is in space 
technology. The installation in the Space Shuttle Columbia of a remote controlled 
manipulator to place satellites into orbit and retrieve them when they fail is just one 
example [9]. Mechanical manipulators have also been used extensively in undersea research, 
probably even more frequently than in space; the latest example being the robotic unit 
used in the discovery of the Titanic [10]. Robot systems could also be used in hospitals to 
help paralytic people or those who must be in bed after surgery. The household robot is 
another dream. Military applications are also appealing. However, until all the control 
problems are overcome, the domain of applications of robot systems will remain limited. 
This will be discussed later, but, for now, a typical structure of a robot system is presented. 

D. STRUCTURE OF A TYPICAL ROBOT SYSTEM 

As illustrated in Figure 1.3, a robot system is, functionally, made up of four interactive 
parts [11]. These different parts are: 
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desired 

evolution control actaul 




Figure 1 .3: A general structure of an adaptive robot 



6 



1. The mechanical manipulator 



The mechanical manipulator, itself, or the plant, is an open chain of rigid links of the 
type previously shown in Figure 1.1. This is the part of the machine designed to perform a 
specific task. Each link is powered by an actuator which physically moves the link in 
accordance with some prescribed control law. The joints are usually equipped with sensors 
to allow for the relative positions of the adjacent links to be measured. 

2. The environment 

The environment in which the robot operates is the physical universe surrounding the 
mechanical manipulator. It includes not only geometrical considerations, but also the 
physical laws governing this universe, the medium in which the robot is immersed and their 
effects on the movements of the robot. Moreover, the robot payload changes constantly 
either by handling parts of different masses or changing tools and configurations from one 
task to the other. This modifies the mass and inertia of the robot, which in turn affects its 
dynamic behavior. These changes must be taken into account in the formulation of the 
dynamic model of the mechanical manipulator. 

3. The task 

The task to be carried by the manipulator, or the trajectory planner, may differ from 
one application to another. In most cases, however, the ultimate task is driving the end 
effector of the manipulator to a desired position in the workspace. Naturally, this position 
is expressed in cartesian or world coordinates. Theoretically, the task might be 
accomplished in any fashion, as along as the tool reaches the final desired position and 
remains there. More realistically, the robot must meet certain requirements in performing 
its task [12]: 

1. The motion must be as fast as possible, otherwise, the use of robots would not have 
been efficient. 

2. No overshoot of the final position is allowed to prevent damages to the environment. 
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3. The mechanical manipulator must be able to avoid obstacles that may be present in its 
workspace. 

4. The motion must be smooth, in order to avoid the increased wear on the mechanism 
and the resonances caused by vibrations due to rough and jerky motions. 

Therefore, stating the initial and final conditions alone is not a sufficient task 
description. In most cases, it is necessary for each link to follow a prescribed trajectory in 
terms of position, velocity, and acceleration at each instant of time. The desired trajectory 
can be made by the combination of any smooth functions joining the initial and final 
positions and satisfying all the constraints. Cubic functions are among the most commonly 
used trajectories since they are easy to generate [13]. 

4. The controller 

The controller generates the control signals that excite the corresponding actuators to 
produce the torques necessary to maintain a prescribed motion of the arm along the desired 
trajectory. The control strategy is determined according to both the control task and the 
mechanical manipulator equations of motion. It is, however, usually derived on the basis of 
a trajectory expressed in joint coordinates. Therefore, a transformation from world 
coordinates to link coordinates must be performed before the signals sent by the trajectory 
planner can be used by the controller. 

In practice, the four functions described above closely interact with each other. When 
in operation, the computer receives, at each instant of time, information concerning the 
robot and information concerning the environment. By using this information in 
conjunction with the control law, it causes the manipulator to move toward the correct 
execution of the task assigned to it. 

This thesis addresses the design of control systems in order for the manipulator to 
adapt to a changing environment, as described by functions 2 and 4 above. The next 
section defines in more detail the mechanical manipulator control problem and addresses 
some of the difficulties. 



8 



E. THE ROBOT CONTROL PROBLEM 

Robotics, while bringing together many well established fields of engineering, is 
relatively a new science in itself. It still suffers from many unsettled points. Controlling the 
robot system to perform in certain way is one of its most challenging problems due to the 
fact that these systems are highly nonlinear. A formal statement of this problem is not, 
however, as difficult as trying to find a satisfying solution for it. In general terms, the robot 
control problem can be formulated as follows: 

Given a desired trajectory generated by the trajectory planner and a mathematical model 
of the mechanical manipulator and its interactions with the environment, find the control 
algorithm which sends torque commands to the actuators in order to cause the desired 
motion to be realized. 

One may now recognize that the robot control problem as stated here is basically a 
stability problem, along the given trajectory. 

Mechanical manipulators may be modeled precisely enough, since their behavior is 
described by the known laws of mechanics [14]. This knowledge should be used in the 
control synthesis as extensively as possible. As stated earlier, this problem is extremely 
difficult because the robot systems are inherently characterized by nonlinear dynamics that 
include nonlinear couplings between the variables corresponding to different motions. 
Furthermore, the dynamic parameters of the manipulator vary with position of the joint 
variables, which themselves vary in time and with respect to each other. These difficulties 
make the implementation of real time dynamic control computationally impractical in 
today's computers. Therefore, one of the intriguing questions arising in the solution of the 
control task is to what extent one should take into account real robot dynamics in control 
synthesis. 

Current industrial practices, in order to take advantage of the well-established linear 
systems and control theory, model the manipulator as a chain of constant-parameter, 
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uncoupled linear subsystems. These design procedures, which may be referred to as the 
servomechanism control methodologies, while yielding satisfactory performances at low 
speeds, have proven to be inefficient for faster and more accurate robotics applications [15]. 

Recently, more researchers have turned to adaptive control in an attempt to be able to 
take advantage of the full robot dynamics and to overcome the limitations of the actually 
available practices [16,17]. These new approaches may be referred to as adaptive control 
strategies for mechanical manipulators. However, no completely acceptable answer has 
been given to the question of how to use the knowledge of the robot dynamics to synthesize 
such control that would be simple enough to implement in practice and to guarantee 
satisfactory system behavior. Several other problems remain [18], such as: 

1. The lack of adequate sensors for the acquisition and pre-processing of information 
received from the environment, particularly visual information. 

2. The state of development of overall theory is not yet fully developed; and, 

3. The slowness of the computations involved. 

The first problem is more of a technological problem than theoretical one and will not be 
dealt with in this context. The last two problems are inherently related to the robot control 
problem and will constitute an important part of this research. 

F. LITERATURE REVIEW 

This section presents some of the most representative solutions to the mechanical 
manipulator control problem as of today. The main difficulty, however, in trying to review 
the literature, is that different approaches have been developed for different classes, types, 
configurations, and purposes. We will primarily consider the approaches to the control 
synthesis for industrial manipulators. We will also restrict ourselves to dynamic control 
which takes into account dynamic effects of the robotics system. Control strategies both in 
terms of open loop and closed loop control have been examined [19]. 
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1. Open loop control systems 



In this case, the trajectory is preplanned or prerecorded and the input torques do not 
depend on link position and/or velocity measurements. The performances are defined in 
terms of desired cost minimization criteria. 

Along these lines, Kahn [20] has considered the time optimal control problem, 
Whitney [21] has studied the minimum energy trajectories, and Young [22] has been 
interested in minimizing a quadratic function in acceleration. 

Due to the highly nonlinear model of the manipulator, numerical solutions only can 
be obtained and stored in memory. In general, this yields a control which is optimal 
provided the system is not affected by unexpected disturbances. Also the open loop 
approach leads to schemes which are very sensitive to parameter variations. Disturbance 
rejection and position tracking can only be achieved through accurate mechanical design. 
The performances of the systems are limited by the capabilities of the actuators and by the 
vibrations induced in the mechanism by the excitation of high frequency structural modes. 
The on-line implementation of such control laws is very involved and might demand a 
rather complex multiprocessor. 

2. Clospd loon control systems 

These feedback control strategies are derived either through well known classical 
servomechanism procedures, or through more recent adaptive control techniques for their 
ability to account for parameter uncertainties. 

a. The servomechanism approaches 

Kahn and Roth [23] have proposed an approximated optimal law which, for a 
particular robot, has resulted in response times and trajectories reasonably close to the 
optimal solutions. For more complex manipulators, however, this solution might be 
unacceptable. Furthermore, the controller proposed in [23] is based on a bang bang 
approach, often unacceptable due to continuous chattering of the joint actuator's signals. 
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Vukobratovic and Stokic [24] addressed the more general problem of designing a 
controller which yields desired tracking while, at the same time, minimizes an appropriate 
cost criterion. 

Because of analytical and computational complexity, approaches by optimal 
synthesis have been developed for positional control problems only. To solve the problem of 
tracking a prescribed trajectory, Popov and co-authors [25] introduced an alternative 
approach which consists of calculating, off line, the nominal trajectory by some optimal or 
suboptimal procedures and then following the obtained path. 

The design of control systems based on the exact nonlinear model of the 

manipulator, in general, yields algorithms not suitable to real time implementations. For 

this reason, controllers based on linearized models in the neighborhood of operating 
conditions have been introduced. This, however, guarantees the stability of the linearized 
model only. Instabilities might occur in the actual system due to nonlinearities in the 

mechanism, coupling between different joints, or parameter variations. In order to 

overcome this major difficulty, several additional compensation schemes have been 
proposed [26,27] at the expenses of added complexity. 

In a different context, Paul [28] has investigated the so called inverse problem 
technique (also named the computed torque by Bejczy [29]). This approach uses the desired 
position, the desired velocity, and the desired acceleration to compute the driving torques. 
The main drawback of this scheme is that the computation of the complete nonlinear 
dynamic model is required. Simplifications have been obtained by Paul [30], Bejczy [31], 
Raibert and Horn[32] by omitting some of the terms in the model. These simplifications, 
while reducing the computational complexity, are still not enough for real time 
implementation of any control strategy based on this technique. 

Vukobratovic and Stokic [33] have recognized that the forces (moments) acting on 
the robot joints can be directly measured and used to synthesize a feedback law that 
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compensates for the coupling in the robotics manipulator and relieves the controller from 
on line computation of these complex terms. Other attempts to include force feedback 
control account for the work of Hewit and Burdess [34]' who introduced force transducers in 
the joints of the manipulator. Although the computation time is shorter, their scheme is 
still too complex for real time implementation. Wu and Paul [35] implemented an analog 
force feedback loop on a single joint manipulator which avoided the computational 
difficulty. Luh, Fisher, and Paul [36] have analyzed the effects of linear independent joint 
torques control. The stability of the overall system, however, has not been discussed in any 
of these papers. There were other attempts to use force feedback, not only at the executive 
control level but to include assembling tasks, such as in the resolved motion control 
introduced by Whitney [37]; the resolved acceleration control by Luh, Walker, and Paul 
[38]; and the resolved force control by Wu and Paul [39]. 

The simplest and most widely used control method today is based on decoupling 
and joint control. Yuan [40] tried to dynamically decouple a manipulator by linear control. 
An effective analysis of a constrained linear control may be found in Golla, Garg, and 
Hughes [41]. Freund [42] attempted the decoupling by nonlinear control involving full state 
feedback which guarantees stability in the absence of external disturbances. Young [43] 
developed a variable structure controller for manipulators, 
b. Adaptive techniques 

In addition to the computational complexity, the servomechanism approaches 
cannot always satisfy the stability conditions even if designed to be robust with respect to 
parametric and state disturbances. Adaptive control methodologies aim at overcoming 
these difficulties. 

Within the adaptive control theory, two fundamental approaches exist in the 
literature [44]. The first is the Learning Model Adaptive Control (LMAC), in which an 
improved model of the plant is obtained by on line parameter estimation techniques, and is 
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then used in the feedback control. A general structure of this approach is shown in Figure 
1.4. The estimated model and the controller may be either linear or nonlinear depending on 
the estimation technique used. The well known Self Tuning Regulator method belongs to 
this class. The second approach in adaptive control theory is the Model Reference Adaptive 
Control (MRAC). The controller is adjusted so that the dynamics of the closed loop system 
matches that of a preselected model. A general structure of this methodology is given in 
Figure 1.5. In general, the reference model is chosen to be a stable, linear, time— invariant, 
decoupled system. The controller may be either linear or nonlinear. It is also possible to 
design adaptive schemes which combine both techniques. 

Many different structures of self tuning regulators are available in the literature, 
and they differ in parameter estimation technique and control algorithms [45]. Koivo and 
Guo [46] examined the feasibility of least squares techniques to robotics applications. Their 
approach is configuration dependent and does not account for nonlinearities in the system. 
It is based on estimation of the linearized dynamics and does not take advantage of any a 
priori knowledge of the system that might be available to the designer. Elliot, Depkovich, 
and Drapper [47] gave an extension of this method to the nonlinear case taking advantage 
of the fact that, in spite of their nonlinear nature, the parameters in the dynamic equations 
of a robot system appear linearly. This method showed better tracking ability, but did not 
solve the computation complexity. Cristi, Das, and Loh [48] exploited this idea of linear 
parameterization of the dynamic equations to give one of the first attempts to formulate an 
adaptive version of the computed torque technique. Their scheme consisted of an on line 
loop to estimate the payload and an off line loop to estimate the other parameters of the 
manipulator. The good feature of this approach is that it guarantees global stability of the 
system. In similar fashion, Craig, Hsu and Sastry [49] proposed another adaptive computed 
torque controller version, based on linear parameterization of the dynamic equations, and 
have established global convergence for their scheme. This method assumes no a priori 
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Figure 1 .4: A Learning Model Adaptive 
Control (LMAC) structure. 
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Figure 1 .5: A Model Reference Adaptive 
Control (MRAC) structure. 
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knowledge about the system and requires acceleration measurements, which makes it 
numerically complex. Middeleton and Goodwin [50] gave an extension of this method based 
on position and velocity measurements only. This method is still complex for real-time 
implementation. Lee and Chung [51] exploited the self— tuning regulator structure by 
introducing the adaptation at the level of linearized perturbation equations in the vicinity 
of a nominal joint trajectory. Their approach uses the recursive Newton-Euler equations 
for feedforward computation of nominal control and a recursive least squares, one step 
ahead control for feedback corrections about the nominal trajectory. The number of 
computations involved is reduced. However, this method can only compensate for small 
deviations. An attempt to speed up this method by avoiding the use of the Newton— Euler 
recursion was performed by Yukobratovic and Kircanski [52]. Their scheme is based on 
local adjustable controllers at each joint. The controller consists of a nominally tuned 
feedforward PID structure, and a feedback corrective portion to account for parameter 
variations. This method relaxes the computational burden by introducing a computer for 
each link, instead of one main computer for the whole robot system. 

There are four basic approaches to the design of Model Reference Adaptive 
Control Systems [53]: 

1. Local parametric optimization theory 

2. Lyapunov functions 

3. Variable structure systems 

4. Hvperstability and positivity concepts 

Within the local parametric optimization techniques, Dubowsky and DesForges 
[54] used the steepest descent method to develop one of the first contributions in adaptive 
control for robot manipulators. This method is computationally less burdensome and has 
good noise rejection properties. However, the steepest descent algorithm, while it can yield 
better adaptation speed, calls for many simplifications and may negatively influence the 
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overall stability of the manipulator. The input signal may also become excessively large 
due to the fact that only the output error is minimized. The discrete time version of this 
method [55], as well as the multivariable case [56], have also been developed by Dubowsky. 
The later approach was tested on an industrial robot and showed the significance of 
adaptive control in high speed tracking operations. 

Takegaki and Arimoto [57] have considered the applicability of model reference 
adaptive control theory in robotics, using the Lyapunov function approach. Their scheme 
included a nonadaptive gravity compensation loop. This resulted in simple adaptation and 
control laws, thus making it suitable for real time implementation. However, how the 
gravity compensation loop affected the tracking quality could not be shown. 

Young [58] combined the variable structure theory with the model— following 
approach and investigated their use in robot positioning problems. This approach also uses 
less computation. It, however, uses the hierarchical control methodology of Utkin [59], 
which may not be valid if only asymptotic convergence can be reached. In addition, the 
control signals are discontinuous and the high frequency components may become 
unacceptable. This method suffers also from the fact that there are no design procedures for 
tuning the controller parameters. Slotine and Sastry tried [60] to remove some of these 
difficulties by using the concept of time varying sliding surfaces in the state space. They, 
however, trade off accuracy against chattering by approximating the obtained 
discontinuous control law by a continuous one. 

Horowitz and Tomizuka [61] presented one of the first attempts to apply 
hyperstability theory to robotics. Their algorithm has been later implemented on a three 
degree of freedom manipulator by Anex and Hubbard [62] after been slightly modified to 
compensate for gravity. The advantage of this method is that the adaptation mechanism is 
derived from the condition of overall system stability. Its main problem is the fact that the 
dynamical effects are estimated without using any a priori knowledge about the system 
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dynamics. In practice, many of the robot parameters are known and it would be convenient 
to estimate only the unknown ones. The application of hyperstability theory to robotics 
models has been fully developed by Balestrino, De Maria, and Sciavicco [63]. Their strategy 
offers better transient behavior compared to that with self— tuning regulators and it 
guarantees stability of the entire system. The main drawback of this method is the 
possibility of excessive control signals and its high numerical complexity. 

As a conclusion, let us note that adaptive methods for manipulation robot are still 
in their early stages of development. It is, therefore, very difficult to produce an exhaustive 
survey of these methods as new design ideas continue to appear in the literature. So far we 
have summarized some of the approaches available, far from a complete treatment of the 
problem. 

G. THESIS OUTLINE 

The remaining of this research will be centered around the development of adaptive 
control strategies applied to robotics systems. Fundamental to the problem of dynamic 
control, the derivation of the dynamic equations of motion will be addressed in Chapter 
Two. Different approaches to obtaining these equations, as well as their computational 
complexity will be discussed. 

In Chapter Three, a new adaptive control law which will combine properties from both 
the STR technique in [51] and the hyperstability principal in [63] will be presented. This 
methodology has the advantage of assuring global stability and aims at overcoming many 
of the limitations of the previously studied schemes. It makes use of a nominal dynamics 
feedforward compensation loop, but, unlike Ref. [51] , the stabilizing feedback loop does not 
call for any simplifying assumptions. A rapprochement between this method and the 
AMFC will be established, but, unlike Ref.[63], it does not require excessive actuation. 
Most of all, our approach yields better performance and is numerically very efficient. 
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II. MECHANICAL MANIPULATOR DYNAMICS 



A. INTRODUCTION 

In order to design a controller of an articulated mechanical system, it is necessary to 
have a mathematical model the system. This model expresses the relationships among 
different components of the robotics system and the interactions between the mechanical 
manipulator as a whole and the physical universe surrounding it. It is described in terms of 
characteristic variables which are specific to the system, such as degrees of freedom, 
lengths, masses, inertias, positions, forces, and torques. 

The number and nature of the parameters used in each model depend on the application 
and the accuracy required. The designer is constantly faced with the challenge of 
developing models that adequately represent the dynamics of the system, and that are 
computationally convenient for computer implementations. 

Because of the high speeds required in any future robotics application, dynamic 
phenomena, such as frictional, inertial, centrifugal, and coupling forces should be taken into 
consideration for the chosen model to be representative of the actual mechanical 
manipulator behavior. 

An efficient mathematical model of a robot is essential for both design and control 
purposes. In the design phase, a complete dynamic model is useful for determining loads, 
dimensions, tolerances and actuation. In control applications, the dynamic model is used to 
generate the nominal joint torques as well as to simulate and test control strategies without 
the need of building a prototype (at least in the early stages of the design). 

There are two problems related to the dynamics of a manipulator. In the inverse 
dynamics problem, we are given a trajectory in terms of joint coordinates q(t) and their 
derivatives. q(t) and q(t), and we wish to find the corresponding sets of vector torques r. 
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This formulation is at the basis of the control problem. The second formulation is the 
direct dynamics problem. In this formulation, we wish to calculate the resulting motion of 
the manipulator q(t), q(t) and q(t) for every given set of vector torques r. This is at the 
basis of the simulation of the robotics system. 

There are several possible approaches one can take to derive the dynamics equations of 
an articulated mechanical system. Newton— Euler's equations, Gibbs' functions, 
d'Alembert's formalism, Bond graphs and Lagrange equations are only few of these 
methods. 

Lagrange and Newton— Euler equations are, however, the most frequently used in the 
literature. In this chapter, we will present several alternative formulations of these two 
methods, address their computational performances, and show that one can easily be 
derived from the other. 

B. CLOSED-FORM LAGRANGIAN MANIPULATOR DYNAMICS 

This formulation was first applied to open loop kinematic chains by Kahn [64] from the 
more general linkage problem of Dicker [65], and has served as the standard manipulator 
dynamics for over a decade. We begin this derivation by presenting the notation used 
throughout the development. 

The links of a manipulator are numbered consecutively from 1 to n starting from the 
base to the tip. By convention, the reference frame is numbered as link 0. The joints are 
numbered so that the joint i connects link i— 1 to link i. An orthogonal coordinate system is 
fixed in each link as follows: 

Z: is directed along the axis of joint i+1, 

x- lies along the common normal from z ^ to z-, and 

y - completes the right handed coordinate. 
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The relative position of two adjacent links is completely described by: 



a -, the distance between the origins of coordinate systems j-1 and i measured along x-, 
Sj, the distance between x-_ ^ and x ■ measured along 

a ., the angle between the and z- axes measured in a righthand sense about x-, and 

0-, the angle between the and axes measured in the righthand sense about z-_y 



This notation is summarized in Figure 2.1. If the joint is rotational, the joint variable will 



be 0 •; if translational, the joint variable will be s, The symbol q . will designate the 



variable for joint i whether it is s- or 0.. The vector q = 

l V 



Jl 

q 2 



represents the generalized 



coordinates of the manipulator and completely specifies its position. In the subsequent 
development, lower case and uppercase regular letters will be used indifferently to 
designate scalar quantities, lower case bold letters to designate vectors, and capital bold 
letters to designate matrices. Subscripts refer to the physical location of the variable, 
superscripts to the coordinates frame the variable is expressed in. Either of these is omitted 
when referring to the base coordinates frame. 



The Lagrange equations for a nonconservative system are: 

d dL 8L 
3t 






= T i' 



i=l,2,...,n 



( 2 . 1 ) 



where, 

L = K - P is the Lagrangian function, 

I< is the total kinetic energy of the manipulator, 

P is the total potential energy of the manipulator, 
q . is the generalized coordinate of the manipulator. 

V 

q • is the first time derivative of the generalized coordinate, and, 
t . is the torque applied to the system at joint i. 

V 
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Figure 2. 1 : The standard axes definitions for connected links 



23 



To find the kinetic energy of the physical system, we need to know the velocity of each 



joint. Let p • = 

L 



4 

r 



denote the coordinates of link i in the reference frame of the link; 



p -, the same point p*- with respect to the base coordinates frame; = 

X X X 



R ^ 1 A ^ 1 

x i 



0 0 0 1 



the homogeneous coordinate transformation which relates the displacement of the x^ 1 link 
coordinate frame to the ( i-1 link coordinate frame; and T^- the coordinate 
transformation which relates the i^ 1 coordinate frame to the base coordinate frame. The 
rotational transformation R- and the translational transformat ion A • are given by: 



i 



2—1 

Ri = 



and, 



*r‘- 



The variables p^ and p* are related by: 



cos#- — sin#-coso- sin#sinu- 
i x x xi 

sin#- cos# -coso- — cos#-sina- 

x xx xx 

0 sina- cosa- 

x x 



1 for a revo 1 ute joint 



( 2 , 2 ) 



’ p a - cos# - ’ 

r X X 


r 


p a- sin#- 

r X X 


II 


U J 


t 



(2.3) 



where 



P; = T Pi 



t, = t? i-L/iy 1 

l X l 

A 



(2.4) 



(2.5) 



Assuming rigid body motion, all the points p- will have zero velocity with respect to the 

X 

th x 

i coordinate frame. The velocity of expressed in the base coordinate frame is: 



v ■ = 
x 



x i 



- dt (p P “ 3t (T 2 P f) 



( 2 . 6 ) 
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Since p . = 0, v • can be written as: 



v i 



i 



<9T. 



- q,- 

• i do ■ 1 

■j = i 4 ; ■ 



p? 



where 



rrpir-1 

— - — = QjTp 1 , 

a q . 



and 



<9T- 

t 



T 1 T 2- T H QjTf 1 -^ 1 , for j < i 



The matrix Q^- being: 






0 - 1 / 0 0 
MOO 
0 0 0 // 
0 0 0 0 



, with 



f or j > i 
i=l ,2,. . .,u 

//=1 and //=0 for 
a revo 1 lit e joint 

u=0 and /*=1 for 
a pr i smat i c joi nt 

<9T- 



(2.7) 



( 2 . 8 ) 



(2.9) 



( 2 . 10 ) 



In order to simplify notations, let us define U - = , then equation (2.9) can be 

tJ dq 



V 



written as follows: 



V 



T j_l QjTf 1 , for j < i 



, j'=l,2,...,n 



( 2 . 11 ) 



for j > i 



Using this notation, v . can be expressed as: 



v i 



Y U • q • 
L ij h j 

j = 1 



( 2 , 12 ) 



The matrix is the rate of change of the point p*- on link i relative to the base coordinate 
frame as q^ changes. It represents both the linear and angular velocities of the link. 
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The kinetic energy of an infinitesimal mass dm on link i is found as: 

dk i = 2 ( x i + Vi + T dm 

dk, = j Tr(v j vt)dm 
where a trace of an nxn matrix A is defined as: 



(2.13) 



n 



Tr(A) = l 



a-- 

n 



i= 1 

Substituting equation (2.12) for v., the kinetic energy of the infinitesimal mass becomes: 

i i 



dk,= * T r 



l l U iA(l’i dm Pi 

X=lk=l 



(2.14) 



For the whole link, 



K r/ dk i 



T r I 

K i 



Tr 



l 1 

L A=1 k=l 



(2.15) 



The integral term inside the bracket is known as the pseudo inertia matrix of all the 
points on link i with respect to the proximate joint of link i expressed in the i* 1 link 
coordinates system. 



y.= 



y t =/p> 


dm 

i 








’ j x 2 -d m 


j x jV j dm 


ix-z .dm 

’ i i 


I* 


• dm ' 

i 


fxjyjdm 


l^ dm 


hi z 2 dm 


\ y 


■ dm 

i 


ix-z dm 

’ i i 




jz 2 -dm 

J X 


\ z 


■ dm 

l 


fi-dm 

' l 


\y i dm 


i z -dm 

1 i 


|d: 


m 



(2.16) 



The inertia tensor I‘ v of link i about its center of mass in the i^ 1 coordinate frame is 



defined as: 



[< fy^l 


-xx 


[ uv [L k J 


u v_ 



dm 



where the indices u,v,k indicate principal axes of the i coordinate frame and 6 is the 
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2— D Kronecker delta. The pseudo inertia matrix J* can then be expressed in terms of the 
inertia tensor I* as: 

r i , t i , r i 



J* = 

i 



yy 


ZZ J * 


\ l 


m-x 


2 


x .y 


xz 


i e 




i 2 — i* 






i* 


xx yy zz 


1 1 


m-y 


xy 


2 


yz 


l* c 






l l -fl 1 -I* 




l l 


I* 


xx yy zz 


mz 


xz 


yz 


2 


i e 


i x c . 


m-y„ 
l y c • 


m .2 
l c- 


m i 



(2.17) 



where c = 

i 



1 






is the center of mass vector of link i from the f l link coordinate frame 



■th 



and expressed in the i link coordinates system. 

The total kinetic energy K of the manipulator arm can be expressed as: 

n 

K = I K i 

2=1 



K =ll 1 I [TrtUftjjuX^ 



i=l A=l/^1 



(2.18) 



Note that the terms J‘- are dependent on the mass distribution of link i and not on their 
position or rate of motion. Hence, the need to be computed only once for evaluating the 
kinetic energy of the manipulator. 

The potential energy P • of the link i is: 

p i = - m P i = ~ m £ 



Tp 1 

l 



(2.19) 



where g = [ g x g^ g, 0 ] is the gravity row vector expressed in the base coordinate system. 
The total potential energy of the manipulator then becomes: 
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( 2 . 20 ) 



T-p* 

VI 



n n 

? =I p i = ~I m £ 

2=1 2=1 

and the Lagrangian function L then becomes: 

i i [ Tr ( u iA j i u ,i)^ + i m 4 T * p i 



( 2 . 21 ) 



2=1A=1A^=1 2=1 

Performing the differentiation in the Lagrange equations and rearranging, we obtain the 

necessary generalized torque r ? . for joint i actuator to drive the i^ 1 link of the manipulator. 

T = d dL_ _ dL_ 

1 **4, 

n A n A A n 

\= l 1 ^xAx^xW I l l ^ v xkA v xi^rl m xi v x/x 

A =27.-1 A=iA— 1 l—l A =2 

2 = 1 , 2 , . . . , 72 



The above equation can be expressed in a matrix notation as: 

n n n 

rl a #t + I I v a^(+8,-' i=1 ' 2 '-'* 

k= 1 k= l/= 1 

or in more compact form as: 



( 2 . 22 ) 



(2.23) 



r(t) = A 



q(t) 



q(t) + V 



q(t),q(t) 



+ G 



q(t) 



(2.24) 



where 



r(t) = 



T i 
1 

r 
. 0 



is an nxl generalized torque vector applied at joints 



q(t) = 



is an nxl vector of the joint variables of the manipulator, 
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q(t) = 



is an nxl vector of the joint velocity of the manipulator, 




q 



n 



q(t) = 




is an nxl vector of the acceleration of the joint variables, and 



A 




is an nxn inertial acceleration related symmetric matrix whose elements are: 

n 

a ik= I Tr ( U A* J A U Ai) ’ (2.24) 



A=ma x(i,k) 

When i=k, a is related to the acceleration of joint i where r- acts and is known as 
u J i 

effective inertia. When #k, a^is related to the reaction torque induced by the acceleration 
of joint k and acting at joint i (known as coupling inertia), and 

rV i 



V 



q(t).q(t) 




is an nxl velocity re lated vector composed of Coriolis and 
centrifugal forces, where 



and 



n n 

\= l l 'ikWl’ i= 1.2,..„n 

k=ll=l 



(2.25) 



n 

^ ikl~ ^A (2.26) 

A=max(z',/:) 

When k=l. the velocity torques are known as centripetal torques, and when 141 as Coriolis 
torques. Friction torques are also velocity related and can be added to this term as: 

fr =Csgn(q i ) + Vq • 

where 

C is a coulomb friction constant, 
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V is a viscous friction constant and, 
sgn is the sign function, and 



G 






is an nxl gravity loading force vector, where, 




(2.27) 



The dynamic equations of motion as given by equation (2.24) are coupled, nonlinear, 
second order, ordinary differential equations. Notice also that equation (2.24) yields the 
solution of the inverse dynamics problem. For every point (q(t),q(t),q(t)) of a given 
trajectory, it yields the required joint torques vector r. This form allows design of a control 
law that easily compensates for the nonlinear effects. Computationally, however, these 
equations are extremely inefficient as they require: 

1. (ijpi)n^ + (^j^)n^ + (^j^)n 2 (l^)tt multiplications, and, 

2. (— j) 7j 4 + (^-)n 3 + (^)n 2 + (^) n additions. 

for every set point in the trajectory. That is, they are of 0(n 4 ) order of complexity, where 
n is the number of links. 

There are two categories of approaches in trying to implement the closed— form Lagrange 
equations in real time control applications: 



1. Simplifying the dynamics by ignoring the least significant terms and correcting errors 
with some feedback compensation. The simplifying assumptions, however, may not hold 
for all speeds and all ranges of applications. 

2. Precomputing terms in the equations and using a gain scheduling approach. 



C. RECURSIVE LAGRANGIAN MANIPULATOR DYNAMICS 

The main reason for the inefficiency of the Uicker/Khan formulation is due to the fact 
that these equations are closed— form expressions and most of the terms are reevaluated 
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many times. To reduce the complexity of these equations, we need to reexamine the above 
derivation and recast it into a recursive form which is computationally more efficient [66]. 
The generalized driving torques given in equation (2.22) can be expressed as: 



n 



r I • *=i A-.» 



A =i 



(2.28) 



where T^ is defined as: 



A A A 

f X - l U A k%+ l l V XkM < 2 ' 29 > 

k= 1 k= 1 b= 1 

The advantage of the above substitution is that equation (2.29) is never used in the 
computation. More efficient recurrence relations for the velocity T^ and acceleration 



are easilv derived bv straightforward differentiation: 

t a = t a-i t a A_1 

*A = ViTa' 1 + T A-l t ) _1 

t A “ Vl 1 *"' + T A-1 U AA^A 

T A = t A-l T )“ 1+t A-l t A" i+t A-l U AA^A +T A-l ,; ’AA4A +T A-l U AA i iA 



T A _ T A-1 T A 1 + 2T A-1 U AA^A + T A-1 U AAA^A + T A-1 U AA^A 



a-i 



(2.30) 



(2.31) 



(2.32) 



Computing the driving torques as given by equation (2.28), and using the recursive 

relations in (2.30), (2.31). and (2.32), results in an 0 (ri 2 ) order of complexity, requiring: 

1 1 

1. 106 jy/T + 620 7 )U — 512 multiplications, and, 

2. 82n“ + 514n — 384 additions. 

The reduction in complexity comes from the fact that to calculate coriolis and 
centrifugal forces, we only need to calculate instead of all the matrices 

Further computations can be saved by noting that: 



Uv, = U,T 



Xi 



ii A 



(2.33) 



Therefore. the generalized torques equation (2.28) can be written: 
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n 



or 



where 



and 



\= l [WM - m A sU t/ r A p A 
A =i 

n n 

r r Tr < u » M^xl-Wu l m x T iA 



A =i 



A =i 



r.-TKiy^-gUjft 



n 



r \ = Sw» 

A =i 

n 

D = T*J*T* + y T* , 1 tJ + 1 JjT| 

i tit L 2 + 1 A A A 
A=i+ 1 

D = J+ + T* ,D l1 

X XX 2+1 Z+l 



C i= 

A = z 

C. = m-p. + tLc., , 

i i v x z+l z+l 

These recursive relations can be computed as follows: 



For z = 1,2,... n 

1. compute T- by equation (2.30) 

2. compute T-by equation (2.31) 

3. compute T • by equation (2.32) 

4. if z = n, continue. Otherwise, set z'=z+l and return to 1. 

5. compute D by equation (2.36) 

l 

6. compute C ? by equation (2.37) 

7. compute r • by equation (2.35) 

l 

S. if z = 1, stop. Otherwise, set z=z-l, and return to 5. 



(2.34) 



(2.35) 



(2.36) 



(2.37) 
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The advantage of this formulation is that its complexity is now O(rc). It requires: 

1. 830?i — 592 multiplications, and, 

2. 675 n — 464 additions. 

Any other reduction in computational complexity can only be obtained by reducing the 
size of the coefficients in' the above complexity polynomials [67]. This can be achieved 
through reformulating the Lagrangian dynamics in terms of 3x3 rotational matrices rather 
than 4x4 rotation— translation matrices. Because 3x3 matrix multiplications require 27 
multiplications while 4x4 matrix multiplications require 64, we get a greater than 50% 
reduction in the coefficients of the computational cost terms. 

The matrix T- relating the orientation of the coordinate system i— 1 and i is now 

V 

j 1 

reduced to a 3x3 rotation matrix R. . A point on link i, expressed in homogeneous 



th x 

coordinates with respect to the i coordinates frame, is now represented as p^ = 



x i 

% 



with p ? the same point with respect to the base coordinate frame. The following sets of 
vectors are also needed throughout this derivation : 



a). the joint i coordinate origin expressed in the j coordinates frame, and, 

7 ih 

Cj, the link i center of mass expressed in the j coordinates frame. 



The quantities, p > o^ and p‘> are related by: 



p • — o • 4 T p ■ 

y i i i y i 



The velocity of p- is then given by: 

t 



p.= v • = 6 • + T -p • 

y i i i i y i 



and the kinetic energy for a particle on link i is as given earlier: 

dk. = ^Tr(v.vT)dm 

dk j = ^Tr(6 1 oT + 2T 1 p^T t p‘pt f TT)dm 



(2.38) 

(2.39) 

(2.40) 

(2.41) 
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(2.42) 



Integrating over all particles in link i, the total kinetic energy LL of link i is given by: 

K • = i Tr(m o oX + 2T noX + T jjTT) 

i 2 v i i i ill ill ' 

where 

nj = f p*dm (2.43) 

Therefore, the total kinetic energy for all the links is given by: 

n 

K = \l Tr(»jM[ + 2T A n)oJ + iyjfa) (2.44) 

A=1 

and the potential energy becomes: 

n 

r = l 'W) < 2 «) 

A=1 

Proceeding as in the homogeneous coordinate formulation, we can write the Lagrange 
equations for each link as: 

a m.' ;it. _ no 

(2.46) 



d dK dK , d? 

T r ati: + — 



aq, dq t aq ; 

which take into account the fact that the potential energy depends on position only. 
Performing the differentiation and rearranging terms yields: 

.Aar , TT tA" V ^ Al 



n 



T , = l f Tr(m A^I + ^ n ^ i I +U Ai I ‘A ii I +U A, J A T I-” l A* T APA 

X % 



X=i 



Recall that 



o x = o • + T o 



A “ ”i 
do 






A T T i 

= U i, 0 A 



aq, 

U Ai= U « T l 

Substituting the above relations into equation (2.47), we obtain: 
n n 

r r Tr<U i! I ("A®A®I + 0 W^ +T >^l +T AJ^)HU tf J m A T>) 

A=i X—i 



(2.47) 

(2.48) 

(2.49) 

(2.50) 

(2.51) 
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or 



r j =Tr(U..D j )- S U ! ,C 1 



(2.52) 



where 



C i ~ 1 m A T A p A “ m p\ + T !'+1 C 2 



i + 1 2+ 1 



(2.53) 



A = i 



has the same structure as equation (2.37), except for the difference in dimensionality. For 
the term D we can write the recurrence: 

i 



D i - l < m A°i 6 I + °1"A TT I + + T 1 J A T I> 

A =i 

n 

D r I < rr; : +1 or’+o;: +1 )(m A 6]; + nhTpvr' +1 T' + %^]; + j^p) + n;oT + j'TT 

A — 2+ 1 



D — T 2 . . D • , + o 2 . , e- , n 2 o? + J ? TX 

2 2 + 12+1 2 + 12+12 2 2 2 

where 



(2.54) 



n 

e ,-= I <»A # I + ^P 

A =i 

e i = e i+l + ”>,*1 + ( 2 '55) 

The term also has the same recursive expression as defined earlier in equation (2.32), 
though presently referring to a 3x3 rotation matrix. The 6^ term is given by: 

6 a = 6 a-i- t a 6 ) -1 < 2 - 56 > 

This formulation decreases the complexity polynomial of the recursive Lagrangian 
formulation with 4x4 matrices by more than 50%. It requires: 

1. 412n - 277 multiplications, and 

2. 320n — 201 additions. 

This translates into 2.195 multiplications and 1,719 additions, for a six degrees of freedom 
manipulator, which is well within the capacity of today's microprocessors. 
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D. RECURSIVE NEWTON-EULER’S MANIPULATOR DYNAMICS 



Another possible formulation of the equations of motion is based on the Newton— Euler 
approach. While the Lagrangian dynamics were reworked with some effort into an efficient 
recursive form, this method, naturally, yields a set of recursive equations which can be 
applied to the links sequentially reducing the computational burden to its minimum 
possible. 

In this derivation, each link is considered as a free body accelerating in space and 
obeying Newton's second law of dynamics for linear movement and Euler's equation for 
angular rotation. 



Using the same notation as previously defined, the vectors o ■ (joint's i coordinate origin 

2 —\ 

expressed in the base coordinate frame), o ■ (joint's i coordinate origin expressed in the 
(i- 1)^ coordinate frame), and o-_^ (joint's (i- 1) coordinate origin expressed in the base 
coordinate frame) are related by: 

°j = vi + °r‘ < 2 - s7 > 

If v._j and w._ 1 are, respectively, the linear and angular velocity of the coordinate 
system with reference to the base coordinates, then the velocity of 



coordinate system (x-.y-,z^) with reference to the base coordinates is: 

dof 1 . . 

v ,=3t— +vi x °r + vi 



Thus the acceleration is given by: 



d2 °i 1 i-i 

Van — f *i-i xo i + 2 w i -i x at 



j ir- 1 

d °2 . , t-l\ , • 

bw M (w M xo. )+v ._ 1 



(2.58) 



(2.59) 



A i — 1 

i—l 

in which 2w- .Xj: is the Coriolis acceleration, w- ,(w- ,xo- ) is the centrifugal 

e—l dt ’ i— l v 2—1 i ' & 

acceleration and x denotes external product of vectors. 

The angular velocity of the system ( x-,y ,z •) with reference to the base w- and its angular 

III l 

velocity with reference to the ( 2 — 1)^ coordinate frame w ? ^ * are related by: 
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and 



But, 



w i = w ^ i + * r ‘ 



• • , • i — 1 

w, = w • , + w • 
i j— 1 i 



Hence (2.61) becomes: 



dw M 

• t-1 aw 2 , T- 1 

w i = ar- + w <-i xw i 



dw j 

= Vi + a — + w i-i xw I 



Since an angular motion of link i is about the z ^ axis, then, 



2-1 

w = 
2 



dw 

a - 



2 - ,o- if link 2 is rotational 

1 — Hi 

0 if link i is translat ional 

Thus, 

^ ^ 

2 /<?• if link 2 is rotational 

l- Hi 

0 if link 2 is translat ional 

Combine (2.60). (2.63), (2.64), and (2.65) to yield: 

w ( + 2 . /q. if link 2 is rotational 

1 — 7 ? — 1*1 






w . = 
l 



■l 1 i-i H i 



w f 
l-l 



if link i is translat ional 



and 


( W. , 
2-1 


w - = 


i 


U. , 

i-l 



(2.60) 



(2.61) 



(2.62) 



(2.63) 



(2.64) 



(2.65) 



( 2 . 66 ) 



; 1 1 11 111 (2.67) 

if link i is translatinal 

Returning to equations (2.58) and (2.59), we note that if link i is translational in 

coordinates {z- r y- 1 .z- 1 ), it travels in the direction z- v with joint velocity q , 

relative to link i-1. If it is rotational in coordinates (^_py-_pZ-_j), it has an angular 

velocitv w 2 t Thus, 

2 



do 



2—1 



3T 



V 1 7 / 

w- o- if link 2 is rotational 

2 2 



z i-l*i 



if link 2 is translat ional 



(2.68) 
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(2.69) 



&o l 1 

l 

at* : 



f dw z ’ _/ 

- — xo- ^x(w* ^ xo 1 ■ 1) if link i is rotational 

X X V X X ' 



3T 



i f link i is translational 



z • i o * 

i — 1 H i 

Combine equations (2.60), (2.63), and (2.68) to yield: 

if link i is rotational 



v i 



w i x °i ; + v i~i 



i-l 



z- /q ■ + w xo +v- , if link t is translat ional 

l — l* l l l l — l 

fV 0 !"' + -,•*(», -oj-') + \-l if link i is 



v i 



z, ,q ; + w 



i-l 



(2.70) 



(2.71) 



rotational 

i -iu T "i xc T' + 2 w . x ( 2 ;-/«i> .. ... , . 

■ _ i it 1 1 ilk % i s 

+ w -x ( w • x o • ) + v- , translational 

L X x l l ' %—i 

Equations (2.66), (2.67), (2.70), and (2.71) describe the recursive relations of velocities 
and accelerations between link i—l and i. To derive the dynamic motion of the mechanical 



manipulator from the above kinematics information, each link is considered as a free body 
accelerating in space and obeying Newton's second law of dynamics: 

d ( m i v c .) 

F = — — r — - — m - a , for linear movement, (2.72) 

l 

and Euler's equation: 

d (I • w,) 

N- = — ^ = I -w • + w-x(I.w-), for angular rotation (2.73) 

where, 

F- is the total external vector force exerted on link i at the center of mass c , 

i i 

dc. 

v c = gy is the linear velocity of the center of mass c ■ of link i, 



a„ = is the linear acceleration of the center of mass c of link i , 
c- dt i 

l 

N . is the total external vector moment exerted on link i at the center of mass c and 

i i 

I is the inertia matrix of link i about its center of mass c- with reference to the base 
i x 

coordinates system. 

The other variables are as defined earlier. 
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The quantities w and w • can be computed from (2.66) and (2.67), while v and a can 
* * c i c x 

be derived from (2.58) and (2.59) as follows: 

Replace o*. 1 by c ) , o - by c-, and note that since both (x.,y.,z.) and c*' are fixed on link i, 

l l v X X X %' X 

do* -1 d2 0 ; -1 

~sr~ -°- 

Consequently, the equations describing v c and a are: 

z z 



\ = w i-i x °r 1 + v i-i 


(2.74) 


v = vi»r* + + *i-i 

X 


(2.75) 


The total external force F • and moment N • are those exerted on link 

l X 


i by gravity and 


neighboring links, that is, 




F • = f • — f • . , 

X X z+1 


(2.76) 


N > = ”r n .+i + (°r c 1 H-(°i+i- c iH+i 

N i = n ; " D i+i + (°, - c ;> xF ; - 0 ! _1 xf i+i 


(2.77) 


where 





f ? - is the force exerted on link z by link i — 1, and 

n- is the moment exerted on link i bv link z — 1. 
z 

Since c ; — o ? _j = o* + c 1 -, equations (2.76) and (2.77) can be expressed in recursive 
relations as 

f i= f , + l + F i < 2 - 78 > 

"i = “i+i + °r’ xf i+i + <°r‘ + c i )xF i + n i (2 - 79) 

According to the convention for establishing link coordinate systems for a mechanical 
manipulator, the motion of link z may only be either a rotation in the coordinate system 
(x ? -_^.y ? -_^.z-_p about z-_j axis, or a translation relative to the coordinate system 
{x- v y _ ,.z _ ,) along z-_ r Therefore, if the joint z is rotational, the input torque r- at 
that joint is the sum of the projection of n ? onto the z-_j axis and the viscous damping 
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moment in that coordinate system. If, however, the joint i is translational, the input force 
r- at that joint is the sum of the projection of f onto the z- 1 axis and the viscous 

X XXI 

damping force in that coordinate system. That is, 

njC z- i + bo- if link i is rotational 

i i — l i^i 

ft z. f + b -o ■ if link i is translat ional 

i i-i i^i 



H 



(2.80) 



where b- is the viscous damping coefficient for joint i. 

In summary, the complete set of equations of motion for the mechanical manipulator 
with n joints and n+1 links consists of equations (2.66), (2.67), (2.70) through (2.75) and 
(2.78) through (2. SO) for i=l,2,...,n. Unfortunately, because these equations are referenced 
to the base coordinate systems, the inertia matrix I- is dependent on the changing 
orientation of link i, which complicates the computation. A more efficient technique for 
computing the joint input forces and torques is to have each link's dynamics referenced to 

j 1 

its own link coordinates [68]. This may easily be achieved using the rotation matrix R- 

c 

defined earlier and noting that since each coordinate system is orthogonal, then: 



(Rjr 1 ) -1 = ( rH) t = 



(2.81) 



Instead of computing w., w-, v-, a £ , F-, N-, f., n., and r-, compute R w^, R w-, R v-, R a , 

i i 

17 XX X 

R F -. R N •. R f, R n-, and R r - Hence the complete set of equations of motion becomes: 

l l l l l 

A-l 



R w = 

i 



R ,_/( R w ,_r + z n^-) if link i is 
11 11 u 1 rotational 



(2.82) 






if link i is 
translational 



'Rj_ ; [R ? i w- j+z 0 q-+{R l ■Axzr.'q) if link i is 
11 1 1 u 1 1 1 u 1 rotational 



R w = 



(2.83) 



R l JR 1 1 w ■ ) 

i-P i' 



if link i is 
trans 1 at i onal 
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if link i is 
rotational 



(2.84) 



R V 



(R^w • ) x ( R 1 o 



i i - 1 



+ (R w .)x [ (R w.) x(R o • ) ] 

R !-,<y i + Ri :Vi>+< R S> x(R<0 "> 

+ 2{R l vf i )x(R l i _ 1 z Q q i ) + (R l v? ^ 
x [ (R ! w-)x(R ! oJ ~ J )] 

R\= (R J w.)x(R t cj0 + (RV.)x[(R ? w.)x(R l C ;.)] + R’v. 



if link i is 
transl at ional 

i. 



R 2 F • = m R V 



R‘N,. = (R'l^XR’w;) + (R‘w i )x[(R ! I i U,.)(R ! w i )) 



l l ' 



R‘ f i= R i+i(R i+1, i + i) + R ‘Fi 



RV=R- +1 [R i+1 n i+1 



l 



+(R 2+1 o) _1 )x(R J+1 f.. 



i + 1 



)] 



+ (R‘oj 1 +R i c^)x(R'F i ) + (R f N i ) 



7 i 



i ? 

(R n ) T (R _ ,z n ) + b -q- if link i is rotational 

l (f J l if 



(2.85) 

( 2 . 86 ) 

(2.87) 

( 2 . 88 ) 

(2.89) 

(2.90) 



(R*f-) T (R! .it,,) + b-q- if link i is translat ional 

if 1 J \J if i 

This formulation gives a 60% reduction in computation over the recursive Lagrangian 
formulation. It requires: 

1. 1507? — 48 multiplications, and 

2. 131 n — 48 additions. 



E. CONCLUSIONS 

In this chapter alternative formulations for deriving the equations of motion of serial 
link manipulators have been described. The emphasis has been put on real time 
computational complexity in terms of required mathematical operations per trajectory set 
point. One should not, however, be misled by the fact that in the above development, the 
recursive Newton— Euler equations are almost three times more efficient than the recursive 
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Lagrangian equations. The discrepancy between the two formulations is due to the 
difference in the angular velocity vector representation used by each method. In the 
Newton— Euler formulation, the angular velocity is adequately represented with a 3x1 
vector w-, whereas in the Lagrangian formulation, it is redundantly represented with a 3x3 
matrix LLy A factor of three in the relative efficiency of the two formulations is therefore 
to be expected. The redundancy in the angular velocity representation is manifested in the 
rotational kinetic energy expression. For the 3x3 representation, the rotational kinetic 
energy is as derived in equation (2.15): 




1 

2 



Tr 



( l l v i^k> 



A=1 



k=l 



whereas for the 3x1 representation, the rotational kinetic energy is: 



K i = 2 w i I i' w i 



(2.91) 



Using this new representation of the rotational kinetic energy, the complete generalized 
force expression r • in the Lagrangian equation (2.22) changes to: 



n 



v I A A “ " 

T i = Z m A(S +t \)-T- + I T A W A + W A X ^A W A^T 



*» A1 



A =i 






dq ,• 



(2.92) 



or 



v r 

i L | f A'^ 

\=i 



A , fVl 

dq i X dqJ 



(2.93) 



where f^ represents the net force in Newton's equation, and represents the net torque in 
Euler's equation. 

Therefore, contrary to what might have appeared earlier, there is no difference in the 
computational complexity of dynamics formulations derived from the Newton— Euler 
equations or the Lagrange equations. The recursive Newton— Euler equations are no more 
efficient than the recursive Lagrangian equations as long as the same representation of 
angular velocity is used. Moreover, the Newton— Euler equations would become as 
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inefficient a s the original Uicker/Kahn equations if they were expressed in closed form. 

Consequently, the emphasis on computational complexity or on advanced control 
strategies synthesis should rest on the structure of the computation rather than on the 
derivation from Lagrange versus Newton-Euler equations. 

In addition, the designer will probably need both structures of the dynamic equations 
and more than one method of obtaining these equations throughout the different phases of 
the design process. He will need: 

1. A closed form expression of the manipulator dynamics in the early stages of the design 
process in order to be able to synthesize adequate control laws, 

2. More than one method of deriving the system dynamics equations in the computer 
simulation phase in order to be able to compare the solution obtained by different 
methods and place greater confidence on the simulation program, and 

3. A recursive form expression of the manipulator dynamics when implementing the 
chosen control law in real time. 
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III. ADAPTIVE CONTROL FOR MECHANICAL MANIPULATORS 



A. INTRODUCTION 

The problem of controlling articulated mechanical systems such as manipulators using 
conventional control methods is very difficult when high speed and high accuracy 
operations are desired. The difficulty arises from the fact that such linkages are 
characterized by highly nonlinear and coupled ordinary differential equations. Closed form 
analytical solutions to these differential equations are not available. Instead, they must be 
solved by numerical integration on a digital computer, which, on the other hand, imposes a 
serious limitation on the number of calculations that can be performed in real time. 

The problem becomes even more difficult when the plant parameters are not precisely 
known and vary in time, as in most robotics applications. Furthermore, a joint angles to 
end point coordinates matrix transformation are usually required in such systems, which 
increases the burden on the computing machine. 

To maintain good performance over a wide range of motions and payloads, researchers 
have turned to adaptive control methods for their ability to adjust to parameters 
uncertainties and load disturbances. Unfortunately, most of these methodologies are not 
computationally efficient. As we have indicated in Chapter 1, two approaches to adaptive 
control theory can be found in the literature. 

In the Model Reference Adaptive Control scheme, the manipulator dynamic model is 
not directly used in the design so that the on line solution of differential equations is not 
required in the implementation. The manipulator is controlled by adjusting position and 
velocity feedback gains to follow a prescribed reference model. 

In the Learning Model Adaptive Control method, a model of the plant is obtained by on 
line parameter estimation techniques. This estimated model is then used in the feedback 
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control. For reasons of computation efficiency, most of the techniques in this category are 
based on linearized models of the manipulator which constrain the range of validity and of 
acceptable performances. Our main focus in this chapter will be to develop an algorithm for 
the Model Reference Adaptive Control of mechanical manipulators. 

In the next section,’ we will formulate the dynamic equations for mechanical 
manipulators in state space. This representation is more suitable for analyzing the 
performances of such methodologies. We will also give a detailed derivation of the state 
space equations of a two link arm model for illustration. 

Section three will review some leading adaptive control methods. These techniques will 
be simulated on the two link model using IBM/DSL [69]. We will compare their 
performances and point out some of their advantages as well as some of their limitations. 
The aim is to show where the need for better adaptive control strategies is felt. 

In section four, a novel adaptive control law which guarantees global stability and yields 
better performances is synthesized. A rapprochement between this method and the Model 
Reference Adaptive Control methodologies is also established. 

Section five summarizes the main results obtained in this study and outlines some areas 
for future research. 



B. MECHANICAL MANIPULATOR DYNAMICS IN TERMS OF STATE VARIABLES 
The standard inverse dynamics equations describing a mechanical manipulator are given 
in equation (2.24) of Chapter 2 and are reproduced here for convenience. 



r(t) = A 



q(t) 



q(t) + V 



q(t),q(t) 



+ G 



q(t) 



These equations can be rewritten as: 



q(t) = A ' 1 



q(t) 



q(t),q(t) 



-G 



q(t) 



(3.1) 



(3.2) 



r(t)-V 

In order to be able to gain better analytical insight to the control system design, it is 
convenient to rewrite equation (3.2) in terms of state variables. 
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Naturally, if one has no particular reason for other choices, the state variables and the 



input control vectors are, respectively, defined as: 

'q(t) 



x p (t) = 



,q(t) 



u p (t) = r(t) 

where, q(t), q(t) and r(t) are as defined earlier in Chapter 2. 



In addition, the vectors V 



q(t),q(t) 



and G 



q(t) 



can be expressed as: 



q(t),q(t) 



G 



q(t) 



= Vi 



= G, 



q(t),q(t) 



q(t) 



q(t) 



q(t) 



Therefore, equation (3.2) can be rewritten as: 



Xp(t) = Ap(Xp(t),t)X p (t) + Bp(Xp(t),t)U p (t) 



w'here 



Ap(Xp(t)’t) = 



V t} = 3t (X P (t)) = 

0 

A" 1 (q(t))Gi(q(t)) 



q(t) 

[q(t)J 

In 



-A-i(q(t))V 1 (q(t),q(t)). 

0 



B p (X p (t),t) = 



A_1 ( q(t))_ 



Here I n is the identity matrix of order n. 



The matrices Vj 



q(t),q(t) 



and Gi 



q(t) 



are found as follows: 



Notice that the expression of V 
rewritten as: 



q(t),q(t) given in equation (2.22) of Chapter 2 



q(t),q(t) 



= V 2 



q(t) 



v 3 



q(t) 



where V 3 



q(t) 



is given by: 



(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

can be 

(3.11) 
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V 3 



q(t) 



and V 3 



q(t) 



Q 1 tt)q i (t) 

q^q^O 
q,(t)q 00 
q (t)q(t) 

• n 1 

q (Oq (t) 



is a matrix containing all the remaining terms of V 



q(t),q(t) 



V 3 



q(t) 



can be expressed as: 



with 



Therefore, 



V 3 



q(t) 



= v 4 



q(t) 



q i (t)q i (t) 

q,(t)q 2 (t) 

q.(t)q (t) 

• 1 n 

q (t)q(t) 

• n 1 

q (t)q (t) 

n n 



q,(t)o.. 0 

0 q (t) 0 ... 0 



q(0 



v 4 



q(t) 



0 0 . . 



q (t) •- 

• D 

0 0 ... 



q,(t) 

0 

q (t) 



( 3 . 12 ) 



( 3 . 12 ) 



To obtain Gi 



q(t) 



, let 



Vi 



q(t),q(t) 



= V 2 



q(t) 



v 4 



q(t) 



q(t) 



l|G[q(t)]|| 
II q(t) II 



q(t) 



( 3 . 13 ) 



( 3 . 14 ) 
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(3.15) 



g] 



q(t) 



Hence, 



Gi 



q(t) 



= H 



G(q(t)] - Gi[q(t)] 
l|G[q(t)] - Gi[q(t)]|| 

l|G[q(t)]|| 



g? 



q(t) 



(3.16) 



where ||.|| is the Euclidean norm of a vector and H 
[70] defined as: 



H 



= I — 2.u.u T 



q(t) II 

is the Householder transformation 

(3.17) 



C. SIMULATION STUDY OF SOME ADAPTIVE CONTROL METHODS 

In the remainder of this chapter, a two revolute joints arm model is considered. This 
model is shown in Figure 3.1. The equations of motions of this mechanical system are 
derived in detailed in the Appendix. Based on this model, a new adaptive algorithm, as 
well as some of the strategies presented in Chapter 1, are studied here in more detail. 
Considered are the inverse dynamics control technique, the variable structure control 
approach, the model following strategy and the perturbation control theory. The basic idea 
behind all these methodologies is to synthesize a control input r which will force the robot 
to follow the output of a reference model. The behavior we are concerned with here is the 
tracking in real time of desired trajectories. The reference model can be either a stable 
linear time invariant decoupled system as in the Adaptive Linear Model Following Control, 
or a combination of models such as in the Variable Structure Control, or a nonlinear 
nominal model of the plant as in the inverse dynamics control and the perturbation theory. 

1. Inverse Dynamics Technique 



In this technique, the control input Up(t) is chosen as: 



U p «=A 



q(t) 



{qO(t)+K v e(t)+K p e(t)}+V 



q(t),q(t) 



+G 



q(t) 



(3.18) 
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Figure 3.1: A two link Mechanical Manipulator 
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where A 



q(t) 



, V 



q(t),q(t) 



and G 



q(t) 



are some estimates of A 



q(t) 



, V 



q(t),q(t) 



and G 



q(t) 



, respectively. 



This control law consist of two basic loops: 



1. A feedforward component: 




= A 



q(t) 



qd(t) + V 



q(t),q(t) 



+ G 



q(t) 



This component is based on a dynamic model of the manipulator. It compensates for 
the interaction forces among various joints. 

2. A feedback component: 



u pb (t) = A 



q(t) 



K v c(t) + K p e(t 



This component is based on position and derivative feedback. It computes the necessary 
correction torques to compensate for any deviations from the desired trajectory. 



Among the attractive features of this method is the fact that, in principle, it turns a 
nonlinear, coupled mechanical system into a linear, decoupled, and stable system. This can 
be seen by substituting the above torques expression of equation (3.18) into equation (3.1) 
to obtain: 



A q + V + G = A {qd + K y e + K p e } + V + G (3.19) 

where the arguments have been omitted for convenience. 

If A = A. V = V and G = G, equation (3.19) reduces to: 



<i(0 



(*K0 + K v e(0 + K p e(oj = 0 



(3.20) 



Since A 



q(t) 



is always nonsingular, K p and K v can be appropriately chosen so that the 



position error vector e(t) approaches zero asymptotically. This control strategy is 
simulated on the two arm model presented earlier. The desired trajectories are chosen as: 



qd(t) = 270t 2 — 180t 3 
qd(t) = 45 + 270t 2 - 180t3 
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(3.21) 

(3.22) 



These trajectories are shown in Figures 3.2 and 3.3. The gain matrices and K y are 
chosen as: 



K = 
P 



K = 
v 



100 0 

0 100 
"20 0 1 



(3.23) 

0 20 . 

In this case, where we assume complete and exact knowledge of the manipulator dynamics, 
we obtain perfect tracking of the desired trajectories. This can be seen in Figures 3.4 and 

j j 

3.5, where the actual trajectories (q^ and q 2 ) and the desired trajectories (qj and q 2 ) are 
virtually indistinguishable with the tracking errors (e^ and e 0 ) shown in Figure 3.G. Also 
the magnitudes of the torques applied at the joint actuators are bounded with reasonable 
values as shown in Figure 3.7. 

The main drawback of this method is inherent in its assumption that one can 



accurately compute the counterparts of A 



q(t) 



, V 



q(t),q(t) 



and G 



this is not always the case in robotics applications. When A 



q(t) 



q(t) 

, V 



. Unfortunately, 



q(t),q(t) 



and 



q(t) 


are not equal to A 


q(t) 


, v 


q(t),q(t) 


and G 


q(t) 



, the quality of the tracking 
degrades and the system may even become unstable. Simulation results with 10% error in 
the load are reported in Figures 3.8 through 3.11. In Figures 3.8 and 3.9, we can see that 
the actual trajectories (q^ and q 2 ) diverge from the desired trajectories (q^ and q^). The 
position error is in the order of 6° in the first link and of 13° in the second link (Figure 
3.10). The applied torques (Figure 3.11) are of reasonable values. 

One may try to overcome this limitation by combining the above scheme with an on 



line identification algorithm to compute A 



q(t) 



, V 



q(t),q(t) 



, and G 



q(t) 



[71]. This, 



however, is computationally demanding since the computed torque, in itself, requires a 
number O (n 4 ) of calculations, n being the number of links in the manipulator. 
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Position (rad.) 




Time (sec.) 



Figure 3.2: The first link test trajectory 
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Position (rad.) 




Figure 3.3: The second link test trajectory 
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Position 




Figure 3.4: The first link desired (q, (t) ) 
trajectories under the Inverse 
(Perfect Modeling) 



and actual (q,(t)) 
Dynamics Controller 
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Position (rad.) 




Time (sec.) 



Figure 3.5: 



The second link desired (q^(t)) and actual (q a (t)) 
trajectories under the Inverse Dynamics Controller 
(Perfect Modeling) 
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Position error (rad.) 



x 1 0 -7 




Time (sec.) 



Figure 3.6: The joint one (e, (t)) and the joint 
tracking errors under the Inverse 
Controller (Perfect Modeling) 



two (e^(t)) 
Dynamics 
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Torques (N/m) 




Figure 3.7: The total torques applied to joint one (£, (t) ) and 
to joint two (t) ) under the Inverse Dynamics 
Controller (Perfect Modeling) 
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Position (rad.) 



1.6 




Figure 3.8: The first link desired (q^(t) ) and actual (q,(t)) 
trajectories under the Inverse Dynamics Controller 
(Approximate Modeling) 
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Position (rad.) 




Figure 3.9: The second link desired (q^(t)) and actual (q^(t) ) 
trajectories under the Inverse Dynamics Controller 
(Approximate Modeling) 
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Position error 



0.1 



0.05 

0 



£ -0.05 
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e , (-0 
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Time (sec.) 



Figure 3.10: The joint one (e, (t) ) and the joint two (e <> (t)) 

tracking errors under the Inverse Dynamics 
Controller (Approximate Modeling) 
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Torques (N/m) 




Figure 3.11: The total torques applied to joint one (*&, ( t ) ) and 
to joint two (^(t)) under the Inverse Dynamics 
Controller (Approximate Modeling) 
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2. Variable Structure System 



A different approach which yields a simple and robust control can be obtained by 
using a variable structure strategy. The theory of variable structure control has mostly 
been developed in the Soviet Union over the last 25 years and has found applications in 
many industrial processes [72]. The fundamental idea behind the theory of variable 
structure is to allow the controller to switch between different strategies, according to 
appropriate functions of the trajectory error. 

A variable structure control is of the form: 



problem consists of choosing the functions u* •, u* •, and the switching hyperplane matrix 

jj t IJ L 

C = diag(c •) such that the sliding mode occurs on the switching hyperplanes, the tracking 
error has an acceptable transient response and it goes to zero asymptotically as t -* oo. This 
methodology is simulated on the two link arm model. The desired trajectories are the same 
as before. The switching planes are chosen as: 



In addition, to ensure the existence of the sliding modes, the control law is chosen as: 



In the absence of a procedural method to selecting the parameters ai, a possible choice 
which will facilitate the calculations is: 



for i = 1, 2,..., n, 




S -(e^) = cf- + e-, c ■ > 0, is the i component of the switching hypersurfaces. The design 



■th 



S j — 0.5e ^ + e ^ 
= 0.4e^ + e^ 




o 1 = 0 , a 2 = 0 , ft 3 = 0 , a 4 = 0 , a 5 = 0 and a = 10 
i ii l i i 



a 1 = 0 , a 2 = 0 , o 3 = 0 , a 4 = 0 , a 5 = 0 and cr = 20 
2 2 2 2 2 2 
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The results of this computer simulation are reported in Figures 3.12 through 3.15. In 
Figure 3.12, the first link actual trajectory (q^) tracks the desired trajectory (q^) with 
approximately 7° error because of poor choice of the parameters oi. The actual trajectory of 
the second link, however, shows better tracking as seen in Figure 3.13. Figure 3.14 shows 
the time evolution of the’errors e^ and e 9 in both joints. As expected, Figure 3.15 shows 
considerable chattering in the input signals. These simulation results highlight the fact that 
chattering in the input signals, absence of a procedural method to choosing the control 
parameters oi, and the difficulty of guaranteeing the existence of the sliding modes are the 
main reasons that limit the applicability of this scheme to multivariable control systems. 

3. Adaptive Linear Model Following Control 

This scheme is depicted in Figure 3.16. The reference model is chosen to be a stable, 
linear, time invariant and decoupled system as: 



X (t) = A X (t) + B U (t) (3.25) 

nr ' mm' m nr ' v > 

where the torques U m (t) are selected so that the output X m (t) of the model follows 
precisely the desired trajectories, described by the user. A m and B jn are of the form: 



A 



m 




(3.26) 



with 





(3.27) 



a „• > 0 , a f • > 0 
m Ox ’ mil 

For the purpose of simulation. a m ^ t - = 1-5 and a m u = 2.5, i — 1, 2. The manipulator 



m 1 1 



input U is chosen as: 
P 



U p (v.X p .t) = *(v.X p ,t)X p - K p X p + *(v,X p ,t)U m + K u U m 



(3.2S) 
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Position (rad.) 




Figure 3.12: The first link desired (qf(t)) and actual (q,(t)) 

trajectories under the Variable Structure 
Controller (Perfect Modeling) 
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Figure 3.13: The second link desired (q^(t)) and actual (q^(t)) 

trajectories under the Variable Structure 
Controller (Perfect Modeling) 
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Figure 3.14: The joint one (e,(t)) and the joint two (e^(t)) 

tracking errors under the Variable Structure 
Controller (Perfect Modeling) 



66 
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Figure 3.15: 



The total torques applied to joint one (^(t)) and 
to joint two ( r S^,( t ) ) under the Variable Structure 
Controller (Perfect Modeling) 
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Figure 3.16: A Linear Adaptive Model Following Control 

for mechanical manipulators 
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where K p and K u are feedback constant gain matrices designed for specific nominal values 
of the plant to satisfy the perfect model following conditions given by: 



K p = -A*(q(t)) (A m - Ap ) 
K u = A '(q(t)) B m 



(3.29) 

(3.30) 



with A + (q(t)) being the pseudoinverse of A(q(t)). 

The quantities $ and ^ are generated by the adaptation mechanism to guarantee the 
stability of the overall system. Possible choices are: 

* = c— 



v 

v 



(sgn(X p ))T 



(sgn(U m )) T 



(3.31) 

(3.32) 



with 



c> 



a 



[ Amax (1UU)] 
Ami n ( A' 1 ( q(t))) 

[ Amax(SS" r )] 
Ami n (A* 1 (q(t))) 



where 



R = A*i(q(t))A*(q(t))( A m - A p ) + A-i(q(t))K p 
S = A-i(q(t))A + (q(t))B m -A-i(q(t))K u 



(3.33) 

(3.34) 

(3.35) 

(3.36) 



Simulation results of this technique with K p and K y obtained from equations (3.29) and 
(3.30) at t = 0 as: 



K p = K p(° } = 



-.354 30.4 6.77 1.72' 

30.4 4.91 1.72 .833 
f2.71 .687 



K u = K u (0) = 



.687 .333 



are reported in Figures 3.17 through 3.20. Since the matching matrices K p and K u are not 
adjusted in time, we can see that the actual trajectories (q^ and q 0 ) follow the desired 



69 



Position (rad.) 




Figure 3.17: The first link desired (q^(t)) and actual (q,(t)) 
trajectories under the Adaptive Model Following 
Controller (Matching Matrices not adjusted) 
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Position (rad.) 




Figure 3.18: The second link desired (q^(t)) and actual (q^(t) ) 
trajectories under the Adaptive Model Following 
Controller (Matching Matrices not adjusted) 
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Time (sec.) 



Figure 3.19: The joint one (e, (t) ) and the joint two (e^(t) ) 

tracking errors under the Adaptive Model Following 
Controller (Matching Matrices not adjusted) 
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Torques (N/m 



250 




Figure 3.20: The total torques applied to joint one ( ^ (t) ) and 

to joint two ('Si(t)) under the Adaptive Model 
Following Controller (Matching Marices not adjusted) 
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trajectories (q^ and q 9 ) with an error of 4° in the first link (Figure 3.17) and of 7° in the 
second link (Figure 3. IS). The joints errors (e^ and e 9 ) are shown in Figure 3.19. The 
control signals (Figure 3.20) are chattering due to the high frequency component. 

Using this type of adaptation, on line numerical integration of the dynamics 
equations of motion is avoided. However, the signals that the actuators are required to 
generate (Figure 3.20) are about 10 times larger than in the computed torque (Figure 3.11). 
This is a serious threat to the plant hardware since the forcing signals are discontinuous. 

To be able to reduce the parameters ( and £, and hence, the actuation signals, one 
should calculate and K u that will satisfy the perfect model following conditions of 
equations (3.29) and (3.30) at each instant of time. Simulation results obtained using this 
fact are given in Figures 3.21 through 3.24. In this case, the trajectory of the first link 
(Figure 3.21) as well as the trajectory of the second link (Figure 3.22) show very close 
tracking, the position errors in both the first and the second links are reduced to zero 
(Figure 3.23). The actuation signals (Figure 3.24) are still large. The main disadvantage of 
this choice is, however, the added computational complexity. 

4. Adaptive Perturbation Control 

A block diagram of this scheme is shown in Figure 3.25. This methodology uses an 
available nominal model of the system and the recursive Newton— Euler equations of 
motion to compute nominal control inputs for a given trajectory. These nominal torques 
compensate for all the interaction forces among various joints along the nominal trajectory. 

To compensate for small deviations from the nominal trajectory, a feedback adaptive 
component is introduced. This adaptive control is based on linearizing the manipulator 
dynamics equations in the vicinity of known nominal trajectory set points to obtain the 
associated perturbed state equation: 

e = Ae + Bdr (3.37) 
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Figure 3.21: The first link desired (qf(t)) and actual (q,(t)) 
trajectories under the Adaptive Model Following 
Controller (Matching Matrices adjusted) 
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Position (rad.) 




Figure 3.22: The second link desired (q^(t) ) and actual (q 4/ (t) ) 
trajectories under the Adaptive Model Following 
Controller (Matching Matrices adjusted) 
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Figure 3.23 



The joint one (ej (t) ) and the joint two (e^(t) ) 
tracking errors under the Adaptive Model Following 
Controller (Matching Matrices adjusted) 
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Figure 3.24: The total torques applied to joint one (^(t)) and 

to joint two (^(t)) under the Adaptive Model 
Following Controller (Matching Marices adjusted) 
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ROBOT LINK PARAMETERS DISTURBANCES 




Figure 3.25: Adaptive Perturbation Control 
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where e and e are as defined earlier and dr = U — U , U n being the nominal torque inputs 
as obtained from certain available nominal model of the manipulator. The system 
parameters A and B depend on the instantaneous manipulator position and velocity along 
the nominal trajectory. A recursive least squares parameters identification technique is 
used to estimate the unknown elements in A and B. The obtained parameters are then used 
to formulate a one step optimal controller that will generate the torques dr to compensate 
for the perturbations. When only the feedforward torques are implemented, simulation 
results show that 10% error in the load produces tracking errors in both the first (Figure 
3.2G) and the second (Figure 3.27) links. As seen in Figure 3.28, these tracking errors are in 
the order of 6° and of 9°, respectively. The input torques (Figure 3.29) stay within 
reasonable limits. Figures 3.30 through 3.33 give the results of the same simulations as 
above when both the feedforward and the correcting torques are used. While an improved 
tracking is experienced in both the first (Figure 3.30) and the second (Figure 3.31) links, 
the joint errors (Figure 3.32) are still of the order of 5° and 3°, respectively. The input 
control signals ( shown in Figure 3.33) stay within the same range of values as before. 
These results are, however, expected since this strategy assumes slow variations and small 
deviations about the desired trajectory. 

It is evident from the above discussion that in order to extend the capabilities of 
manipulators and improve their overall dynamic performances, there is a need to 
investigate and develop better adaptive control solutions to current control problems. 

The aim of the next section is to present a novel adaptive control law for mechanical 
manipulators that enjoys global stability and overcomes some of the limitations of the 
previously studied methodologies. This strategy combines properties from both the the Self 
Tuning Regulator in [51] and the Model Reference Adaptive Control in [63] and offers itself 
to microcomputer implementation. This technique serves also to extend the Model 
Reference Adaptive Control method into using a nonlinear reference model. 
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Position (rad.) 




Figure 3.26: The first link desired (q^(t)) and actual (q,(t)) 

trajectories under the Adaptive Perturbation 
Controller (No Adaptation) 
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Figure 3.27: The second link desired (q^(t) ) and actual (q^(t)) 

trajectories under the Adaptive Perturbation 
Controller (No Adaptation) 
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Figure 3.28: 



The joint one (e,(t)) and the joint two (e^(t)) 
tracking errors under the Adaptive Pertubation 
Controller (No Adaptation) 
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Figure 3 . 29 : The total torques applied to joint one (2T ( (t)) 
and to joint two (2*,(t)) under the Adaptive 
Perturbation Controller (No Adaptation) 
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Figure 3.30: The first link desired (q^(t)) and actual (q,(t)) 

trajectories under the Adaptive Perturbation 
Controller (With Adaptation) 
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Figure 3.31: The second link desired (q^(t)) and actual (q^(t) ) 

trajectories under the Adaptive Perturbation 
Controller (With Adaptation) 
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Figure 3.32: The joint one (e, (t) ) and the joint two (e^(t)) 

tracking errors under the Adaptive Pertubation 
Controller (With Adaptation) 



87 



Torques (N/m) 




Figure 3.33: The total torques applied to joint one ('£,( t) ) 
and to joint two ('^(t)) under the Adaptive 
Perturbation Controller (With Adaptation) 
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D. ADAPTIVE NONLINEAR MODEL FOLLOWING CONTROL 

Figure 3.34 illustrates the structure of the proposed adaptive control system. The task of 
the controller is to generate the control signals in order to follow a desired trajectory 
despite the changes in the manipulator's parameters and the errors in the dynamic model. 
The desired trajectory is specified in terms of joint angles q d (t) and their derivatives q d (t) 
and q d (t). The total torques r are obtained through a nominal and a correction loops. 

The nominal loop is justified by the fact that, in practice, a nominal model of the 
manipulator is always available to the designer. The nominal parameters are used to 
construct a recursive inverse dynamics that generates nominal torques r . The inputs q n , 
q , and q j( to the recursion are obtained by adding filtered error signals z(t), z(t), and z(t) 
to the desired signals q d , q d , and q d . The signals z(t) are chosen such that: 

Z(s) = F[s)E{s) (3.38) 

where F[ s) is the transfer function of a linear filter to be determined and Z(s) and E(s) are 
Laplace transform of z(t) and e(t) respectively. The recursive form is chosen to ease the 
computation. The inverse dynamics form is motivated by the fact that, in the case where 
the mechanical manipulator model is known precisely, an exact trajectory following is 
obtained. The adjusting signals z(t) are selected to reach the ideal closed loop dynamics 
given by the error equation: 

e(t) + K y e(t) + K p e(t) = 0 (3.39) 

Since some manipulator parameters such as load and inertia vary in time and some 
others such as friction in the gears and motors backlash are very difficult to determine, 
additional correction feedback torques dr are necessary to account for any deviations from 
the desired trajectory due to these effects. These correcting torques are generated to 
guarantee global stability: 

lime( t) = 0 (3.40) 

t-* oc 
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Figure 3.34: Adaptive Nonlinear Model Following Control 

( Series Configuration ) 
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(3.41) 



In total, the input torques r axe the sum of r Q and dr. 

t = t + dr 
n 

where r Q are the outputs of the recursive algorithm given by: 



r D W = A 






q n M + v 



q n (t),q n (0 



+ G 



q„« 



(3.42) 



and dr are to be determined. From the equations of motion of the manipulator, the total 
torques r can also be expressed in terms of the joint angles q(t) as: 



r(t) = A 



q(t) 



q(t) + V 



q(t),q(t) 



+ G 



q(t) 



(3.43) 



Subtracting equations (3.43) from (3.42) and substituting dr by its expression in equation 
(3.41) yields: 



q(0 



or, since A 



q n (0 -q(t) 



q(t),q n (t),q(t),q n (t) 



q(t) 



= w 

is nonsingular, 

q n -q = -A' 



+ dr 



W + dr 



where the arguments have been omitted for convenience and where, 



W=AV 



with 



q(t),q Q (t),q(t),q n (t) 



+AG 



q(t)>q n (t) 



+AA 



q(t)iq n (t) 



m 



AV = V 



q^ij.q^t) 



- V 



q(t),q(0 



AG = G 




-G 


q(t) 


AA = A 


q n « 


-A 


q(t) 



lim W 



q(t)'q n (t) 5 q(Oiq n (t) 



= o 



as q(t) — q Q (t) and q(t) — q fl (t). 
On the other hand, 



(3.44) 



(3.45) 



(3.46) 

(3.47) 

(3.48) 

(3.49) 

(3.50) 



q n - q = q d + £-'(s J f(s)fl(s))-q 



(3.51) 
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That is, 



q n -q = e + L*i(s 2 i^s)£(s)) (3.52) 

where is the inverse Laplace transform of the given function. Substituting the 

expression of (q fl — q) from equation (3.52) into equation (3.45) we obtain: 



e + L~ l (s 2 F\s)E(s)) =-A'i W + dr (3.53) 

Equation (3.53) shows that the adjusting forward signals are crucial to the stability of the 
system. When these signals are not used, such as in [51], the error equation is unstable. 

Any stable filter F{ s) of the form: 

<tn- 2 S n ~ 2 + a n -3S n ' 3 + .. . + ao 



F(s)= (3.54) 

b n s n + b n - iS n ’ 1 + . .. + bo 

will yield the desired error equation in (3.39). A more interesting choice, however, is: 

n s) = — (Ks + K ) (3.55) 

S 2 v P 

where K y and K p are velocity and position feedback matrices, respectively. The motive 
behind choosing F{ s) as in equation (3.55) lies in the fact that it relaxes the computation 
considerably. This can be seen by evaluating z(t), z(t), and z(t) corresponding to this 
choice: 



z(t) = L' l (s 2 F\s)E(s)) 

z(t) = K y e(t) + K p e(t) (3.56) 

which is actually no more than velocity and position measurements feedback. The 
quantities z(t) and z(t) can then be obtained from the above expression of z(t) by a simple 
and a double integration, respectively. With F{ s) selected as in equation (3.55), the error 
equation in (3.53) becomes: 



e + K e+K e = -A‘i 

v p 



W + dr 



(3.57a) 



Equation (3.57a) is an equivalent error model representation of the proposed adaptive 
control law. It can be partitioned into a linear time invariant system connected with a 
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nonlinear time varying block in the feedback as shown in Figure 3.35. It should be clear 
that K y and Kp are chosen such that the forward transmission function is stable. 

Notice that, if a filter of the type given in equation (3.54) is used instead, the term is 
decomposed as: 

F[s) = 1 - F(s) 



Hence, equation (3.53) becomes: 



with 



e + K y e + Kpe = - A ’ 1 



Wi + dr 



(3.57b) 



Wj = W + L-i(F(s)E{s)) 
which is a more general form of equation (3.57a). 

It remains now to determine the control inputs dr based on the knowledge of some 
upper bounds of: 



w |q(t).q n (t),q(t),q n (t) 
such that || e(t) || < e, for t -> » and for any initial conditions. Here c is a small number. 
Let 



is) = s(t) = L-wmm) 



i)(t) = K y e(t) + K p e(t) 



The quantity tj{ t) is given by: 

fa) = - K 



arj + A' 1 



W + dr 



with o=K' ] K = diag(o-), a >0 for i=l,...,n , since 



"v p oV V' i 



K = diag(k .), k • > 0 and K = diag(k -), k • > 0. 

V V V VI p ov p v pz 

Define dr as: 



dr(t) = f(t) 



fa) 

Mil 



with 



(3.58) 

(3.59) 



(3.60) 



(( t) > Amax(A)|| WTA-i|| 



(3.61) 
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(s 2 I + K v s + K p ) 



e(t) 



A M (W + dr) 



Figure 3.35: Error model representation of the proposed 

adaptive control law 
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FACT: 



Equations (3.60) and (3.61) imply: 



Consequently, 



lim T]{t) = 0 
00 



h‘me(t) — 0 

£-*oo 

lim e(t) = 0 

£-*oo 

Zm(q n -q d ) = 0 

£-*oo 



PROOF: 



The above claim can be proven by choosing a Lyapunov function V(r?) as follows: 

v(*rt = -W(tM t) 

Then 

V(iy) = *7 T (t)v(t) 



v(»/) = - 



rjtaK T}+ WTA-iK 77+drTA-iK 17 



V ' V ' V 

Notice that since 0 = diag(ap, > 0 and K v = diag(k vJ ), k V} - > 0, then 



vv' vi 



oK y = 0 = diag(/T), 0 - > 0. 

Now replace dr(t) by their expressions in (3.60) to obtain: 



V{t]) = — 7 f0T}— W t A' 1 K v 77+ Crf^ A' l K^ 



If we select f(t) so that: 



Amin(A-i)£> || WTA'i| 



then 



That is. 



with 



V(t?) < - rp0T) < - J7 T (t)»?(t) 



V(,)<-2/? m . n V(,) 



(3.62) 

(3.63) 

(3.64) 

(3.65) 



(3.66) 

(3.67) 

(3.68) 



(3.69) 

(3.70) 

(3.71) 

(3.72) 
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0mn = ain W j),f° r >= h-,n. 

by which (3.62), (3.63), (3.64) and (3.65) are satisfied. 

The proof for the case where the filter in (3.54) is used, follows along the same lines. 



QED 



The above control law is simulated on the two link arm model. The following values of 

K , and K are used: 
p’ v 



/ = 



16 

0 



V 



K = 
v 



0 ' 

10 . 

0 

4 

0 

4 



When 10% error in the load is assumed, the tracking quality, when no feedback corrections 
are applied, is poor in both the first (Figure 3.36) and the second (Figure 3.37) links. The 
joint errors are of the order of 6° and 17°, respectively as shown in Figure 3.38. The toque 
signals (Figure 3.39) are similar to the ones obtained from the inverse dynamics law. There 
is no noticeable improvement (or little) in the tracking quality of both the first (Figure 
3.40) and the second (Figure 3.41) links, when the adjusting signals z(t), z(t), and z(t) 
alone are used. This can be seen from the plots of the time evolution of the tracking errors 
in Figure 3.42, and of the input torques in Figure 3.43. However, the quality of the tracking 
is improved drastically in both the first (Figure 3.44) and the second (Figure 3.45) links, 
when both the adjusting signals and the correction torques are implemented. Figure 3.46 
shows that the tracking errors in both links are practically reduced to zero. This is about S° 
better than [51] and 10° better than [63]. As expected, figure 3.47 shows a moderate 
chattering, within a very acceptable range of values, in the control signals. Because 
integrators are used to generate z(t) and z(t), it is of interest to check that the presence of 
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Figure 3.36: 



The first link desired (g?(t)) and actual (q, (t) ) 
trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (No Feedback) 
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Figure 3.37: The second link desired (q A (t)) and actual (q.*,(t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (No Feedback) 
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Figure 3.38: The joint one (e,(t)) and the joint two (e^(t)) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (No Feedback) 
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Figure 3.39: The total torques applied to joint one ( & ( (t) ) and 

to joint two (^(t)) under the Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (No Feedback) 
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Figure 3.40: The first link desired ( q , ( t ) ) and actual (q, (t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (Using Adjusting Signals alone) 
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Figure 3.41: The second link desired (q A (t)) and actual (q^(t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (Using Adjusting Signals alone) 
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Figure 3.42: The joint one (e ( (t)) and the joint two (e^(t)) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (Using Adjusting Signals alone) 
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Figure 3.43: The total torques applied to joint one ( 'S’, (t) ) and 

to joint two ("^(t)) under the Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (Using Adjusting Signals alone) 
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Figure 3.44: The first link desired (q,(t)) and actual (q, (t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric Uncertainties 
only (With Adjusting Signals and Correcting Torques) 
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Figure 3.45: The second link desired (q^(t) ) and actual (q*,(t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric Uncertainties 
only (With Adjusting Signals and Correcting Torques) 
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Figure 3.46: The joint one (e,(t)) and the joint two (ej,(t)) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller assuming Parametric Uncertainties 
only (With Adjusting Signals and Correcting Torques) 
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Figure 3.47: The total torques applied to joint one ( £ (t) ) and 

to joint two (?i(t)) under the Adaptive Model 
Following Controller assuming Parametric Uncertainties 
only (With Adjusting Signals and Correcting Torques) 
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offset measurements in such signals does not affect the tracking quality of the overall 
system. Simulations results with initial conditions of e^ = — 115° are given in Figures 

3.48 through 3.51. There is practically no degradation in the tracking quality of either the 
first (Figure 3.4S) or the second link (Figure 3.49). The tracking errors (Figure 3.50) and 
the input torques (Figure 3.51) are as before. To show the robustness of this methodology 
to parameter disturbances, 50% error in the load is assumed and a term of the form F = 
Csgn(q-) + Vq . is added to the plant as unmodelled friction, where C and V are Coulomb 
friction and viscous friction constants, respectively. The simulation results from this case 
with C = V = 4 are reported in Figures 3.52 through 3. 03. The actual trajectories (q^ and 
q. 7 ) diverge considerably from the desired trajectories (q^ 1 and q 9 ) when no feedback is used 
as seen in Figures 3.52 and 3.53, respectively. The tracking errors are in the order of 32° in 
the first link and of 114° in the second link as shown in Figure 3.54. The required input 
torques are reported in Figure 3.55. A considerable improvement in the tracking quality of 
both links is achieved even when only the adjusting signals alone are used as can be seen in 
Figures 3.56 and 3.57, respectively. However, the joint errors are still of the order of 12° in 
the first link and of 23° in the second link as seen from Figure 3.58. The control signals are 
reported in Figure 3.59. When both the adjusting signals and the feedback correcting 
torques are implemented, the tracking quality is close to perfect as can be seen in Figure 
3.60, for the first link and in Figure 3.61, for the second link. The errors in both links 
(Figure 3.62) are drastically reduced to approximately 0.3° and 0.1°, respectively; while the 
input torques (Figure 3.63) remain at very acceptable range of magnitudes with mild 
chattering. When these same conditions are simulated with the Adaptive Perturbation 
Control Law of [51], the quality of the tracking is poor for both the first (Figure 3.64) and 
the second (Figure 3.65) links. Figure 3.66 shows the time variations of such errors. These 
errors are of an order of magnitude of 15° in the first link and of 12° in the second link. 
Respectively, this is about 50 and 120 times less accurate than the proposed approach (see 
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Position 




Figure 3.48: The first link desired (q,(t)) and actual (q,(t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (With 115° offset) 
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Position (rad.) 




Figure 3.49: The second link desired (q^(t)) and actual (q^(t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (With 115° offset) 
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Figure 3.50: The joint one (e, (t)) and the joint two (e^(t)) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (With 115° offset) 
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Torque (N/m) 




Figure 3.51: The total torques applied to joint one (£ ( (t)) and 

to joint two (^(t)) under the Adaptive Model 
Following Controller assuming Parametric 
Uncertainties only (With 115° offset) 
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Figure 3.52: The first link desired (qf(t)) and actual (q,(t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (No Feedback) 
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Figure 3.53: The second link desired (q A (t)) and actual (q^(t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (No Feedback) 
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Figure 3.54: The joint one (e, (t) ) and the joint two (e.j,(t) ) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (No Feedback) 
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Figure 3.55: The total torques applied to joint one (T/(t)) and 

to joint two (^X(t)) under the Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (No Feedback) 
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Figure 3.56: The first link desired (q,(t)) and actual (q, (t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (With d^ = 0) 
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Figure 3.57: The second link desired (q^(t)) and actual (q^(t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (With d£ = 0) 
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Figure 3.58: The joint one (e, (t) ) and the joint two (e z (t)) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (With d'S' = 0) 
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Figure 3.59: The total torques applied to joint one ("&(t)) and 

to joint two ( ^(t)) under the Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (With d£" = 0) 
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Figure 3.60: The first link desired (qf (t) ) and actual (q,(t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (Complete Scheme) 
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Figure 3.61: The second link desired (q^(t) ) and actual (q^(t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (Complete Scheme) 
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Figure 3.62: The joint one (e, (t) ) and the joint two (e^(t) ) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (Complete Scheme) 
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Figure 3.63: The total torques applied to joint one ('S’, (t) ) and 

to joint two ('^(t)) under the Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (Complete Scheme) 
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Figure 3.64: The first link desired (q*f(t)) and actual (q,(t)) 

trajectories under the Adaptive Perturbation 
Controller (With Unmodeled Disturbances) 



126 



Position (rad.) 




Time (sec.) 



Figure 3.65: The second link desired (q^(t)) and actual (q^(t) ) 

trajectories under the Adaptive Perturbation 
Controller (With Unmodeled Disturbances) 
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Figure 3.66: The joint one (e,(t)) and the joint two (e 4/ (t) ) 

tracking errors under the Adaptive Pertubation 
Controller (With Unmodeled Disturbances) 
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Figure 3.62). Also, as seen in Figure 3.67, the input torques are about 10 times larger than 
the input torques required by our control law (see Figure 3.63). Similarly, when the same 
conditions as above are simulated with the Linear Model Following Control Law of [63], 
The actual trajectories (q^ and q.2) diverge from the desired trajectories (q^ and q^) as 
seen in Figures 3.68 and 3.69. Figure 3.70 shows that the tracking errors reach, 
approximately, 29° in the first link and 35° in the second link. Here, again, the superiority 
of the proposed control methodology is evident. 

A deeper insight into this adaptive control method may be gained by considering its 
relation to the general structure of Model Reference Adaptive Systems and to the Adaptive 
Model Following Controller in particular. 

Referring to Figure 3.34, the same system can be represented in a slightly different, but 
equivalent, arrangement as shown in Figure 3.71. This particular representation highlights 
more clearly the parallel structure of this control law. 

Now consider Figure 3.72 which gives a block diagram representation of the standard 
Model Following Control law [73]. A fundamental difference between Figure 3.71 and 3.72 
is the fact that the flow of signals through the reference model in Figure 3.71 is reversed. 
This results from our formulation of the general manipulator control problem in which we 
assume that in practice a desired trajectory is specified by the user and not the input 
torques to the manipulator. 

The reference model in the standard Adaptive Model Following Control is chosen to be 
asymptotically stable. That is, A m is a Hurwitz matrix. The Newton— Euler recursion we 
are using is also a stable algorithm, in the sense that for every desired trajectory point 
(q d ,q d ) where |q d | < oc and |q d | < 00, it yields a vector torque where | r | < 00. 

The feedforward matrix gain K u , the plant feedback matrix gain K , and the model 
feedback matrix gain K m in Figure 3.72 are chosen such that, for null initial conditions and 
specific plant parameter values, perfect model following exists. That is: 
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Figure 3.67: The total torques applied to joint one (£, (t)) 
and to joint two ( ^j,(t)) under the Adaptive 
Perturbation Controller (With Unmodeled Disturbances) 
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Position (rad.) 




Figure 3.68: The first link desired (q^(t)) and actual (q , ( t) ) 
trajectories under the Adaptive Model Following 
Controller (With Unmodeled Disturbances) 



Position (rad.) 




Figure 3.69: The second link desired (qf (t) ) and actual (q^(t) ) 
trajectories under the Adaptive Model Following 
Controller (With Unmodeled Disturbances) 



132 



Position error 




Figure 3.70: The joint one (e, (t) ) and the joint two (e^(t)) 

tracking errors under the Adaptive Model Following 
Controller (With Unmodeled Disturbances) 
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Figure 3.71 : Adaptive Nonlinear Model Following Control 

( Parallel Configuration ) 
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Figure 3.72: The standard Adaptive Linear Model Following Control 
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/m A (t) = A (t) 

t~*CC ^ 


(3.73) 


Hm B (t) = B m (t) 

t-* 00 y 


(3.74) 


lira e(t) = 0 

/->oo 


(3.75) 


lim e(t) = 0 

t~* CO 


(3.76) 



In the proposed adaptive control methodology of Figure 3.71 , the reference model is an 
available nominal model of the plant itself used as an inverse dynamics. That is, in at least 
a limited range of applications for which it has been calculated, this model adequately 
describes the plant under consideration. When this is the case, perfect model following can 
be achieved without the need for any adjustment, by choosing = 0, K m = 0 and K u = 
I n . This particular choice of Kp, K m , and K u means that in the case where the values of 
the plant parameters are precisely known and do not vary during operation, the adaptation 
mechanism is not needed just as in the standard Adaptive Model Following Control. In 
fact, the control law reduces to an inverse dynamics control. Simulation results of this 
situation are shown in Figures 3.73 through 3.76. Figures 3.73 and 3.74 show that, in the 
ideal case where all the parameters are known, the actual trajectories (q^ and q 9 ) follow 
very closely the desired trajectories (q^ and q 9 ), with no need for adaptation. The time 
evolution of the joint errors (e^ and 62) are converging to zero as we can see from Figure 
3.75. Figure 3.76 show's that the corresponding input torques are within an acceptable 
range of values and of reasonable chattering. 

As for the standard Adaptive Model Following Control, when the perfect model 
following exists, the role of the adaptation is to assure the convergence to this solution 
when the plant parameters are uncertain or vary during operation. This is shown to be the 
case in equations (3.64) through (3.70). This adaptation law can be classified as a signal 
synthesis adaptation since the feedback signals z(t), z(t) and z(t) are used to either reshape 
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Figure 3.73: The first link desired (q,(t)) and actual (q ( (t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller (Perfect Modeling) 
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Figure 3.74: The second link desired (q^(t)) and actual (qa,(t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller (Perfect Modeling) 
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Figure 3.75: The joint one (e,(t)) and the joint two (e^(t)) 

tracking errors under the Nonlinear Adaptive Model 
Following Controller (Perfect Modeling) 
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Figure 3.76: The total torques applied to joint one ( 2", (t) ) and 

to joint two (%,(t)) under the Adaptive Model 
Following Controller (Perfect Modeling) 
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the forward torques r n or to generate additional torques d r which are both acting as input 
signals to the plant when the values of its parameters differ from the nominal ones. 

The fact that the flow of the signals through the' reference model in Figure 3.71 is 
reversed, does not, theoretically, affect the equations governing the overall system. This 
can be seen from the equivalent error model representation of equation (3.57) and Figure 
3.35 which is the same structure as for the standard Adaptive Model Following Control. 

Finally note that this method is very insensitive to the parameters ^(t)- Simulation 
results when 



and when 



i(t) = 



<t) = 



'20 

0 

120 

0 



0 ' 
16 

0 

1 IS 



are reported in Figures 3.77 through 3.80. Figures 3.77 and 3.78 are practically the same as 
Figures 3.60 and 3.61 obtained earlier with a different gain matrix ^(t). The robustness 
property is inherent to the fact that the proposed control law is a switching law as can be 
seen from Equation (3.60). Figures 3.79 and 3.80 show that, when the gains matrix is 
excessively large (8 times the nominal values), the quality of the tracking is degraded. 
However, the joint errors experienced in this case may still be more tolerable than the the 
joints errors experienced with the control law of [51] or the control law of [63]. 

It is also very important to point out that the proposed adaptive control law is 
numerically more efficient than [51] since it does not explicitly estimate the feedback 
model, and more efficient than [63] because it does not require the use of model matching 
matrices. This should be very attractive feature since adaptive control techniques suffer 
from computational complexity in general. 
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Figure 3.77: The first link desired (q,(t)) and actual (q,(t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (With Large Gains) 
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Figure 3.78: The second link desired (q^(t)) and actual (q^(t)) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (With Large Gains) 
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Figure 3.79: The first link desired (qf(t)) and actual (q, ( t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (Very Large Gains) 
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Figure 3.80: The second link desired (q (t) ) and actual (q (t) ) 

trajectories under the Nonlinear Adaptive Model 
Following Controller assuming Parametric and 
Unstructured Disturbances (Very Large Gains) 
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E. CONCLUSIONS AND FUTURE RESEARCH 

A new adaptive control law for mechanical manipulators that maintains uniformly good 
performance over a wide range of motions and payloads has been developed. This control 
strategy has been shown to combine properties from both the Model Reference Adaptive 
Control and the Self Tuning Regulator theory. It has also been shown that this method 
serves to extend the Adaptive Model Following Control Approach into using a nonlinear 
model as a reference. 

The design procedure is simple and systematic resulting in an overall system which is 
globally stable and offers itself to microprocessor implementation. We have also shown that 
this control law is robust with respect to variations of the plant parameters. The 
effectiveness of the approach has been demonstrated on several computer simulations which 
compare its performances against some of the commonly known adaptive control 
techniques. In all cases, the proposed adaptive control strategy has performed better. 

This adaptive control scheme has reduced the chattering in the input torques to a 
"reasonable" value compared to [59] and [63]. We are currently investigating ways to 
eliminate this chattering completely. As has been shown in the previous section, the 
chattering is the result of the correction torques attempting to counterbalance the effect of 
errors in the manipulator parameters. If the manipulator model is precisely known, the 
correction torques are reduced to zero and the controller becomes an inverse dynamics. 
This can be achieved by off line identification tests. Also it has been shown that in practice 
most of the manipulator parameters can be measured or estimated beforehand and only the 
parameters that are load dependent are unknown [74]. Using this fact, simulations with on 
line recursive least squares estimation of the load alone are currently underway. 
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APPENDIX 



DYNAMICS OF THE TWO LINKS PLANAR MECHANICAL MANIPULATOR 



A. INTRODUCTION 

The two revolute joints planar mechanical manipulator shown in Figure A.l is used as 
the basis for our simulations throughout this study. In this Appendix, the dynamic 
equations describing the motion of this physical system are derived using the Lagrangian 
Euler equations of Chapter 2. 



B. NOTATION 

The same notation and conventions as established in Chapter 2 are also employed here. 
In addition, for i = 1, 2 , the following variables are used to denote: 

O' the joint angle, which also serves as the generalized coordinate; 
m- the mass of link i; 

X 

l- the length of link z; 

l ■ the distance from the proximate joint to the center of mass of link i; and 

C l 

I - the moment of inertia of link i about the axis z- 
i x 



C. EQUATIONS OF MOTIONS 

Equation (2.15) of Chapter 2 can be used to derive the kinetic energy K • of link i. This 
equation can also be broken down into a translational and a rotational parts as: 

K ; = 4- m ! v ci V C> + 4- W U> w ; 



(A.l) 



where v^. denotes the linear velocity of the center of mass of link i and v/- the angular 
velocity of link i about z-. 
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Figure A. 1 : A two link Mechanical Manipulator 
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The total kinetic energy K of the manipulator is then: 

K=XK. 

2=1 l 

The total potential energy is found using equation (2.20) of Chapter 2 as: 



(A.2) 



P = “ i?, m i g c i 

where c - is the position of the center of mass of link i and g the gravity row vector. 
From Figure A.l, 

-l cl B x sin(^) 
l cl cos(^) 

0 

- {Ij sin (flj) + l c2 sin(^+^ 2 )} ^ — ^sin(^+^ 2 ) # 2 
cos (0j) + /^cos^T^)} e x + l c2 cos(^+^ 2 ) &2 
0 

0 



(A.3) 



v c 1 ~ 



y c2 ~ 



w ; = 



0 

1 



(A.4) 



(A.5) 



(A .6) 



and 



Therefore, 



w 2~ 



0 

0 

^1 + ^2 



1 2 • 
Kj = 0 j 



(A.7) 



( A.S) 



K 2 — — y~ m 2 {If ^ h^2 cos (^2^ ^1 

9 

1 



+ 



1 



2 2) 2 
2~ "> 2 ~i e 2 



m t 



L 



H — }j— nif) l 2 cos($f+ $ 2 ) {lfCos(0f)-\ — 2~~ cos(^^+^)} Bf6r) 
+ ~ y ~ m 2 ^sin(^ + ^ 2 ) {/^sin(^ 1 ) + -7^— sin^T^)} 0j0 2 
+ ~W~ m 2 l 2^1 + ^ 



(A.9) 
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ij 4 

P = (m^ + 2 mg) g -jj — sm(^) + m^g — sin^ + fl^ (A. 10) 

^ l 2 

where l ^ and have been replaced by — and — , respectively. 

Without loss in generality, we assume 1^=1^= /, and perform the operations in the 
Lagrangian equations. After rearranging, we obtain the following actuators torques: 



r ll 




a ll a 12 




i 

r-H 

i 




V lll V 112 v 121 v 122 


1 ” 
to 

1 




a 21 a 22 
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v 211 v 212 v 221 v 222 



*?. 

6 \ e 2 

h e \ 

*2 



+ 



where 



a, , — 



11 
a 12 = 
a 21 = 
a 22 = 



3 

1 



m 



3 

1 



m 



2 



U 

u 



3 

1 



m 



3 

1 



f + 



2 

1 



f 4- m 9 f cos (4) 



cos(# 9 ) 



m. 



? cos(0 9 ) 



v lll = 0 
V 1 1 2 = ~ m 2 
v 121 = 0 



v 122 “ 



V 211 “ : 
v 212 = 0 
V 221 = 0 
v 222 = 0 



m 



Z 2 sin(# 9 ) 



%2 f sin(# 2 ) 
l sin(0 9 ) 



'2 



(A.ll) 



m 1 cos(9, +Q 0 ) 

Gj = g Z { — — cos(9]) + m,2 (cos(#j) H ^ )} 

G 2 = — ^ — 1^2 g f cos(^ 1 +^ 2 ) 

Equations (A.ll) through (A. 25) constitute what is known as the inverse dynamics form 
of the equations of motion of the two link mechanical manipulator of Figure A.l. These 



(A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 
(A- 17) 
(A.18) 
(A.19) 
(A. 20) 
(A.21) 
(A. 22) 
(A. 23) 

(A. 24) 
(A. 25) 
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equations are used in Chapter Three to evaluate the performance of many commonly 
known adaptive control algorithms as applied to robotic manipulators. The reason for this 
choice is that Equations (A. 11) through (A. 25) are relatively simple enough to keep the 
operations manageable, and yet, representative of the systems under study. 
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